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ABSTRACT 


A  study  was  made  of  design  techniques  for  low-density 
wind  tunnels  operating  with  high  stagnation  temperatures  and 
hypersonic  flow. 

An  approximate  boundary  layer  analysis  was  developed  to 

predict  feasible  ranges  of  operating  conditions  for  low-density 
wind  tunnels. 

Detailed  calculations  were  made  for  three  nozzles  includ¬ 
ing  boundary  layer  corrections  for  isentropic  contours. 

It  became  obvious  that  the  use  of  high  stagnation 

temperature  wind  tunnel  for  investigations  of  rarefied  gas 
dynamics  could  produce  situations  in  which  it  would  not 

always  be  clear  whether  the  experimental  results  were 
associated  with  rarefied  gas  flow  phenomena  or  related  to 
high  temperature  physical  phenomena. 

It  appears  that  unless  an  enormous  amount  of  money  is 
devoted  to  the  construction  of  an  extremely  large  low-density 
wind  tunnel,  it  will  never  be  feasible  to  test  large  models 
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such  that  Knudsen  numbers  in  the  free-aolecular-f low  range 
are  produced. 


This  Technical  Documentary  Report  has  been  reviewed  and 
is  approved. 

ft  9+^  \ — 

ROBERT  L.  COLLIGAH,  JR.  VJ  \J 

Colonel,  USAF 

Deputy  for  Test  and  Sqpport 


REVIEW  AKD  AFPROVAL  CP  ASD  TECHNICAL  DCCUNiLlTARY 
REPORT  ASD-TDB-63-154 


PROJECT  ENGINEER: 


PAUL  A.  CZYSZ  w 


Cnief,  High  Temperature 
Aerodynamics  Division 


Section 


CONCURRED 


WILLIAM  J.  DuBOIS 
Lt  Colonel,  USAF 
Ac  ting  Director  Engineering  Test 
Deputy  for  Test  and  Support 


CONCURRED  IK: 


iical  Director 
Deputy  for  Test  and  Support 


APPROVED  BY: 


ROBERT  L.  COLLIGAN,  JRU 
Colonel,  USAF 

Deputy  for  Teat  and  Support 


TABLE  OF  CONTENTS 


Page 


INTRODUCTION  1 

Objectives  of  Study  1 

Rarefied  Gas  Dynamics  3 

Present  Status  and  Problems  6 

Conclusions  8 

THE  LOW  DENSITY  WIND  TUNNEL  9 

Configuration  9 

Static  Pressure  Requirements  10 

Stagnation  Pressure  Requirements  17 

APPROXIMATE  BOUNDARY  LAYER  ANALYSIS  19 

Introduction  19 

Sumnary  of  Existing  Information  20 

Summary  of  Approximate  Boundary  Layer  24 

Analysis 

Discussion  of  Assumptions  of  Approximate  28 

Analysis 

Results  of  Approximate  Boundary  Layer  35 

Analysis 

Conclusions  39 

REAL-GAS  EFFECTS  IN  LOW-DENSITY  WIND  TUNNELS  42 

Description  of  Real-Gas  Effects  42 

Typical  Results  45 

Conclusions  51 

DETAILED  NOZZLE  CALCULATIONS  53 

Introduction  53 

Nozzle  Design  Procedure  53 

Nozzle  Design  I  58 

Nozzle  Design  II  65 

Nozzle  Design  III  70 

Summary  70 


vi 


SUGGESTED  DESIGN  PROCEDURE  FOR  LOW-DENSITY  78 

WIND  TUNNELS 

COMPONENTS  FOR  LOW-DENSITY  WIND  TUNNELS  82 

Mass  Flow  Rate  Requirements  82 

Requirements  for  Heater  84 

Requirements  for  Precooler  85 

Requirements  for  Exhaust  System  86 

CONCLUSIONS  88 

RECOMMENDATIONS  92 

REFERENCES  121 

APPENDICES 

A.  STATIC  PRESSURE  REQUIREMENTS  FOR  LOW  U4 

DENSITY  WIND  TUNNELS 

B.  APPROXIMATE  BOUNDARY  LAYER  ANALYSIS  117 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1  Drag  Coefficient  on  a  Cylinder  at  M  =  5,92  under 

Rarefied  Gas  Flow  Conditions  95 

2  Schematic  of  Low-Density  Wind  Tunnel  96 

3  Reynolds  Number  Requirements  for  Low-Density 

Wind  Tunnels  97 

4  Variation  of  Pressure  with  Mach  Number,  Reynolds 

Number  and  Knudsen  Number  for  Air,  T0  =  7000°R  98 

5  Variation  of  Static  Pressure  with  Mach  Number, 

Reynolds  Number  and  Knudsen  Number  for  Air, 

Tq  =  1 2 , 000°  R  99 

6  Nozzle  Stagnation  Pressure  Requirements  100 

7  Nomenclature  for  Approximate  Boundary  Layer 

Analysis  101 

8  Boundary  layer  Thickness  Parameter  for  a  Sine- 

Arch  Velocity  Profile  102 

9  Comparison  of  Laminar  Boundary  Layer  Thickness 

for  Sine-Arch  Velocity  Profile  and  Flat  Plate 
Solutions  log 

10  Laminar  Boundary  Layer  Thickness  Constant  for 

Compressible  Flow  Over  a  Flat  Plate  I04 

11  Minimum  Exit  Diameter  Reynolds  Number  (Isentropic 

Contour)  for  (  re  =  1  at  Exit  of  Actual  Nozzle  I03 

12  Dependence  of  Exit  Diameter  Reynolds  Number 
(Isentropic  Contour)  on  Isentropic  Core  Size  at 
Exit  of  Actual  Nozzle 

13  Isentropic  Contour  Exit  Diameter,  Displacement, 

Thickness  and  Actual  Exit  Diameter  for  a  Low- 
Density  Wind  Tunnel  Nozzle  I07 

14  Actual  Exit  Diameter,  Boundary  Layer  Thickness, 
and  Isentropic  Core  Diameter  for  a  Low-Density 

Wind  Tunnel  Nozzle  ^qq 

15  Actual  and  Approximate  Nozzle  Expansion  Processes  1  qq 


viii 


LIST  OF  ILLUSTRATIONS  (CONT'D) 


Figure  Pa9e 

16  Variation  of  Composition  for  the  Expansion  of 
Oxygen  in  a  15°  Hyperbolic  Nozzle,  T0  =  5900°K, 

p0  =  0,965  Atm.  HO 

17  Variation  of  Composition  for  the  Expansion  of 
Oxygen  in  a  15°  Hyperbolic  Nozzle,  T0  =  5900°K, 

p0  =  9.4  Atm.  1H 

18  Variation  of  Composition  for  the  Expansion  of 
Oxygen  in  a  15°  Hyperbolic  Nozzle,  T0  =  5900°K, 

p0  =  82  Atm.  J 

19  Variation  of  Static  Temperature  for  the  Expansion 
of  Oxygen  in  a  15°  Hyperbolic  Nozzle,  TQ  =  5900°K, 

p0  =  82  Atm.  H3 

20  Variation  of  Static  Temperature  for  the  Expansion 
of  Air  in  a  Hypersonic  Nozzle,  T0  =  6000°K,  p0  = 

100  Atm.  127 

21  Constant*  Approximations  for  Frozen  Flow  Regions 

of  a  Hypersonic  Nozzle  Expansion  Process  128 

22  Variation  of  Velocity  for  the  Expansion  of  Air  in 

a  Hypersonic  Nozzle,  T0  =  6000°K,  p0  =  100  Atm.  129 

23  Variation  of  Mach  Number  for  the  Expansion  of  Air 

in  a  Hypersonic  Nozzle,  T0  =  6000°K,  p0  =  100  Atm.  130 

24  Variation  of  Static  Pressure  for  the  Expansion  of 
Air  in  a  Hypersonic  Nozzle,  T0  =  6000°K,  p0  = 

100  Atm.  131 

25  Variation  of  Density  for  the  Expansion  of  Air  in  a 

Hypersonic  Nozzle,  T0  =  6000°K,  p0  =  100  Atm.  132 

26  Recommended  Values  for  the  Boundary  Layer  Thickness 

Constant  C  133 

27  Mass  Flow  Rate  Requirements  for  Low-Density  Wind 

Tunnels  Using  Air  134 

28  Heater  Requirements  for  Low-Density  Wind  Tunnels 

Using  Air  135 

29  Volume  Flow  Rate  Requirements  for  Low-Density  Wind 

Tunnels  Using  Air  136 

lx 


NOMENCLATURE 


a 

a 

A 

A 

A 

A 

B 

C 

C 

C 


C' 


D 

D 

D* 

D 


D' 

se 


D 


e 


g 

h 


o 


k 

kx,k2 


velocity  of  sound 

radius  of  throat  section  of  nozzle  in  centimeters. 
Ref.  (20) 

cross-sectional  area  of  nozzle 

constant  in  boundary  layer  correction  method, 

Eq.  (50) 

concentration  of  atoms,  Eq.  (37) 
concentration  of  molecules,  Eq.  (37) 

constant  in  boundary  layer  correction  method, 

Eq.  (51) 

boundary  layer  thickness  constant  based  on 
PeVeLt/uo  '  Eq*  (18) 

boundary  layer  thickness  constant  for  flat-plate 
flow,  Eq.  (19) 

boundary  layer  thickness  constant  for  predicting 
boundary  layer  growth  in  nozzle.  Fig.  26 

constant  in  boundary  layer  correction  method, 

Eq.  (52) 

cylinder  diameter 

diameter  of  nozzle  cross-sectional  area 
diameter  at  nozzle  throat 

diameter  at  exit  section  of  isentropic  nozzle 
contour 

diameter  of  isentropic  core  at  exit  section  of 
actual  nozzle  contour 

diameter  at  exit  section  of  actual  nozzle  contour 
acceleration  given  to  unit  mass  by  unit  force 
enthalpy 

Boltzman's  constant 

constants  in  boundary  layer  correction  method, 

Eqs.  (43)  and  (44) 


x 


Knuds en  number 

frozen  flow  parameter.  Ref.  (20) 
any  characteristic  dimension 
distance  along  flat  plate  surface 

axial  distance  between  nozzle  throat  and  exit 
sections 

mass  of  molecule 
Mach  number 

pressure 

minimum  heater  power 

radius  of  nozzle  cross-sectional  area 

radius  of  isentropic  core  at  exit  section  of 
actual  nozzle 

dimensionless  nozzle  radius,  Eq.  (45) 
individual  gas  constant 
universal  gas  constant 

Reynolds  number  based  on  some  L  and  a  Pu/u  product 
for  a  single  state 

Reynolds  number,  Eq.  (48),  PeugL  t/uo 
Reynolds  number  per  unit  length,  ou/u 

molecular  speed  ratio 

volume  flow  rate  at  nozzle  exit  density,  Eq.  (59) 

absolute  temperature 

velocity 

dimensionless  velocity,  Eq.  (47) 
velocity 

mean  velocity  for  Maxwellian  velocity  distribution 
mass  flow  rate 

molecular  weight  of  dissociated  gas 


WQ  molecular  weight  of  undissociated  gas 

x  dimensionless  distance  along  nozzle  axis,  Eq.  (49) 

y  distance  normal  to  a  surface 

ym  number  of  molecules  per  Initial  molecule  of  air 

y  number  of  atoms  per  initial  molecule  of  air 

fl 

Z  compressibility  factor,  Eq.  (7) 


Greek: 

a 

Y 

6 

6* 

6+ 

e 


A 

u 

p 


mass  fraction  of  gas  dissociated,  Eq.  (37) 
isentropic  exponent 

boundary  layer  thickness  (99%  velocity) 

boundary  layer  displacement  thickness 

dimensionless  boundary  layer  thickness,  Eq.  (41) 

total  angle  for  divergent  section  of  a  conical 
nozzle 

mean- f ree-pa th 

viscosity 

density 

dimensionless  density,  Eq.  (46) 
dimensionless  speed  of  sound  parameter,  Eq.  (6) 


Subscripts: 

a 

A 

B 

C 

e 


atmospheric  conditions 

nozzle  section  at  which  actual  expansion  process 
deviates  from  equilibrium- flow  process 

nozzle  section  at  which  actual  expansion  process 
becomes  a  frozen-flow  process 

nozzle  section  at  which  a  frozen-flow  process  is 
started  on  equilibrium  flow  curve 

exit  section  or  conditions  for  actual  nozzle  contour 


n,n  +  1  successive  sections  along  nozzle 


xii 


o 

r 

s 

8 

se 


nozzle  stagnation  state 
a  reference  condition 
isentropic  process 
boundary  layer  seam 

exit  section  or  conditions  for  isentropic  nozzle 
contour 

free -stream  conditions  for  a  flat-plate  flow 


Superscript; 

*  throat  section  of  nozzle 


xiii 


INTRODUCTION 


Objectives  of  Study 

This  report  summarizes  some  analytical  results  of  a  study 
of  low-density  wind  tunnels.  Specific  operating  conditions  for 
these  wind  tunnels  are  selected  so  that  experimental  data  can 
be  obtained  for  flows  about  models  corresponding  to  the  various 
regimes  of  rarefied  gas  dynamics.  The  general  considerations 
for  the  design  of  such  wind  tunnels  are  presented  along  with 
specific  designs  for  several  wind  tunnels.  The  study  was 
supported  by  the  Aeronautical  Systems  Division  of  the  United 
States  Air  Force  Systems  Command  at  Wright  Patterson  Air  Force 
Base,  Ohio  and  is  part  of  ASD’s  continuing  policy  to  improve 
and  to  obtain  experimental  facilities  which  will  support  a 
high-level  "in  house"  research  program.  Since  the  future 
plans  of  the  group  interested  in  hypersonic  flow  involves  an 
experimental  program  in  rarefied  gas  dynamics,  this  study  was 
initiated  to  determine  the  limitations  on  the  use  of  low  den¬ 
sity  wind  tunnels  in  the  rarefied-gas- dynamics  region.  For 
those  operating  conditions  such  that  the  use  of  the  wind  tunnel 
is  feasible,  specific  design  information  is  desired  for  plan¬ 
ning  new  facilities  at  ASD. 

The  objectives  of  this  study  as  set  forth  in  the  contract 

were : 

1.  To  determine  the  state  of  the  art  on  the  design 
of  low  density  wind  tunnels  which  could  produce  flow 
conditions  such  that  models  of  reasonable  size  could 
be  tested  in  the  transition  flow  and  free -molecular- 
flow  regimes. 

i manuscript  released  by  the  authors  February  1963  for  publication  as  an 
ASD  Technical  Documentary  Report. 
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2.  To  determine  those  ranges  of  operating  conditions 
where  a  low  density  wind  tunnel  could  be  utilized. 

3.  To  determine  which  phases  of  the  nozzle  design 
procedure  must  be  modified  so  as  to  account  for  the 
conditions  encountered  in  low-density  hypervelocity 
wind  tunnels  and  to  generate  these  new  techniques 
insofar  as  time  permitted. 

4.  To  design  several  nozzles  for  use  with  available 
and  proposed  high  temperature  heaters  at  ASD.  A 
range  of  stagnation  temperatures  between  7000#and 
12000° R  was  selected  for  the  study.  The  nozzle  de¬ 
signs  were  to  be  based  upon  the  use  of  air  or  other 
non  toxic  gases  and  the  nozzle  exit  velocities  were 
to  be  as  high  as  possible  consistent  with  the  stag¬ 
nation  temperature  range.  The  nozzle  designs  were 
to  include  contour  specifications  and  the  pressures, 
temperatures,  flow  rates,  etc.,  associated  with 
these  contours. 

5.  To  survey  the  available  literature  relative  to 
rarefied  gas  dynamics  with  a  special  emphasis  on 
free-molecular- flow  testing  techniques. 

The  results  of  the  literature  survey  were  published  in 

it 

Reference  (1)  .  Some  700  references  were  presented  on  rare¬ 
fied  gas  dynamics  and  topics  associated  with  low-density  wind 
tunnel  design.  The  references  were  cross  indexed  on  a  subject 
basis.  The  literature  survey  and  this  report  constitute  the 
final  report  of  this  study. 

*  Numbers  in  parentheses  refer  to  items  in  the  Bibliography. 
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Rarefied  Gas  Dynamics 


In  order  to  discuss  intelligently  the  design  considerations 
for  low-density  wind  tunnels  it  is  desirable  to  review  briefly 
the  basic  factors  which  control  rarefied  gas  dynamics.  This  is 
best  accomplished  by  reference  to  Fig.  1  which  shows  some  re¬ 
cent  experimental  data  on  the  drag  coefficient  for  a  cylinder. 
These  data  were  obtained  by  Maslach  and  Schaaf  (2)  at  the  Uni¬ 
versity  of  California.  The  drag  coefficient  is  plotted  versus 
the  Knudsen  number  (based  on  cylinder  diameter)  and  the  results 
were  obtained  at  a  Mach  number  of  about  6.  According  to  its 
usual  definition,  the  Knudsen  number,  Kn,  is: 

A 


where  A  is  the  mean- free -path  and  D  is  the  cylinder  diameter. 
It  will  be  seen  that  at  small  values  of  the  Knudsen  number, 
about  0.01,  the  experimental  values  of  the  drag  coefficient 
approach  the  value  corresponding  to  continuum  flow.  However, 
at  a  Knudsen  number  of  0.1,  deviations  from  the  continuum  flow 
value  are  already  about  10  percent.  As  the  Knudsen  number  is 
increased,  the  drag  coefficient  deviates  more  and  more  from  the 
continuum  flow  value  until  at  a  Knudsen  number  of  about  10,  the 
drag  coefficient  is  2.8  and  has  practically  achieved  a  value 
corresponding  to  the  free -molecular^ low  limit. 

Normally  the  first  departures  from  the  results  of  contin¬ 
uum  gas  dynamics  occur  at  a  Knudsen  number  of  0.01.  The  flow 
regime  in  the  Knudsen  number  range  from  0.01  to  0.1  is  called 
the  slip  flow  regime.  In  this  regime  the  first  effects  of 
particle-surface  interaction  are  encountered  in  that  it  is  no 
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longer  possible  to  relate  in  a  simple  manner  flow  properties 
in  the  fluid  near  a  surface  to  those  properties  at  the  sur¬ 
face.  In  the  case  of  continuum  gas  dynamics,  values  of  tem¬ 
perature  (or  temperature  gradient),  velocity,  mass  transfer 
rates,  etc.,  in  a  fluid  at  a  surface  are  equated  to  values  of 
these  quantities  at  the  surface.  In  the  range  of  Knudsen 
numbers  between  0.1  and  10  the  particle-surface  interactions 
become  more  significant  and  this  flow  regime  is  called  the 
transition  flow  regime.  Finally  when  the  Knudsen  number  is 
10  or  greater  the  free-molecular- flow  regime  is  encountered 
and  the  particle-surface  interactions  are  controlling.  The 
two  limiting  values  indicated  in  Fig.  1  for  the  free-molecular- 
flow  regime  correspond  to  the  cases  of  diffuse  and  specular 
reflection  of  particles  from  the  surface.  Table  I  summarizes 
the  flow  regimes  and  the  corresponding  Knudsen  number  ranges 
which  are  involved. 


TABLE  I 

FLOW  REGIMES  OF  RAREFIED  GAS  DYNAMICS 


Regime 
Slip  Flow 
Transition  Flow 
Free-Molecular-Flow 


Knudsen  Number  Range 
0.01  -  0.1 
0.1  -  10 
10 


The  results  shown  in  Fig.  1  also  depend  upon  the  value  of  the 
Mach  number,  M,  defined  as: 

M-I  (2, 

or  equivalently  the  molecular  speed  ratio,  S,  defined  as: 
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(3) 


I2  S0  RI  ' 

The  relationship  between  the  Mach  number  and  the  molecular  speed 
ratio  is: 


where  Y  i-8  the  ratio  of  specific  heats. 

Based  on  the  above  discussion  it  is  apparent  that  in  the  de¬ 
sign  of  low-density  wind  tunnels  it  is  necessary  to  select  tunnel 
operating  conditions  such  that  correct  values  are  produced  for 
those  independent  variables  which  control  the  heat,  mass  and  mo¬ 
mentum  transfer  processes  of  the  three  flow  regimes  in  rarefied 
gas  dynamics.  In  addition  to  the  obvious  need  for  certain 
values  of  the  Knudsen  number  and  the  molecular  speed  ratio,  it  is 
necessary  to  control  all  factors  which  enter  into  the  particle- 
surface  interaction  especially  in  the  transition  and  the  free- 
molecular-flow  regimes.  It  is  easier  to  state  that  those  inde¬ 
pendent  variables  for  the  particle-surface  interaction  be  con¬ 
trolled  than  to  identify  them.  These  independent  variables  in¬ 
clude  the  surface  condition  of  the  model,  the  temperature  of  the 
surface  and  the  temperature  of  the  incoming  stream.  However, 
even  use  of  the  term  temperature  for  the  incoming  stream  is  a 
substitute  for  an  exact  specification,  since  we  are  more  prop¬ 
erly  concerned  with  the  distribution  function  of  the  incoming 
stream  and  with  the  equilibrium  or  non- equilibrium  state  of  this 
stream.  In  the  results  that  follow  some  of  these  variables  will 
be  related  directly  to  the  tunnel  design;  some  of  the  variables 
will  be  introduced  in  a  somewhat  vague  manner  and  some  of  them 
will  be  no  more  than  commented  upon. 
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Present  Status  and  Problems 


At  the  present  time  most  of  the  experimental  work  on  ex¬ 
ternal  aerodynamic  problems  in  rarefied  gas  dynamics  has  been 
obtained  with  either  of  two  types  of  experimental  equipment. 
The  first  piece  of  equipment  is  the  low-density  wind  tunnel. 
The  second  piece  of  equipment  is  the  molecular  beam  apparatus 
for  studying  particle-surface  interaction  phenomena.  Most  of 
the  existing  experimental  data  from  low-density  wind  tunnels 
have  been  obtained  with  relatively  low  tunnel  stagnation  tem¬ 
peratures,  the  order  of  a  few  hundred  degrees  Fahrenheit,  and 
in  the  Mach  number  range  from  0  to  approximately  7.  It  is 
only  within  the  last  two  years  that  some  experimental  drag 
data  for  spheres  were  obtained  in  a  hypersonic  wind  tunnel  (3) 
and  to  the  authors'  knowledge  these  are  the  only  experimental 
data  taken  at  a  high  stagnation  temperature .  Therefore,  there 
is  a  need  for  additional  test  facilities  for  obtaining  aero¬ 
dynamic  data  on  models  of  various  shapes  under  high  stagnation 
temperature  and  high  Mach  number  flow  conditions  in  the  rare¬ 
fied-gas-dynamic  regions.  However,  a  word  of  caution  is  re¬ 
quired  on  the  use  of  such  high  temperature  hypersonic  facili¬ 
ties  In  general,  a  high  temperature  plasma  tunnel  tends  to  be 
a  difficult  piece  of  apparatus  in  vAiich  to  measure  a  large 
number  of  individual  quantities.  Although  the  plasma  tunnel 
serves  an  important  function  in  producing  experimental  data 
which  can  be  subsequently  scaled  for  use  in  the  solution  of 
design  problems  for  re-entry  vehicles,  etc.,  it  is  not  the 
best  piece  of  equipment  from  which  to  get  accurate  and  suffi¬ 
cient  data  to  confirm  or  deny  proposed  theoretical  results. 

For  this  reason  it  is  well  to  keep  in  mind  that  there  will  be 
a  continuing  need  for  a  variety  of  low-density  wind  tunnels 
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capable  of  operation  over  a  wide  range  of  stagnation  tempera¬ 
tures  and  Mach  numbers,  and  using  gases  which  are  not  nec¬ 
essarily  related  to  any  practical  problem,  present  or  future. 
This  requirement  exists  because  it  is  necessary  to  identify 
those  factors  which  are  related  to  rarefied  gas  dynamics  and 
those  factors  which  are  related  to  high- temperature  physical 
phenomena . 

The  available  analytical  and  experimental  information 
on  the  surface-particle  interaction  problem  of  rarefied  gas 
dynamics  is  poor.  A  recent  survey  of  accommodation  coeffic¬ 
ients  by  Harnett  (4)  lends  independent  support  to  this  con¬ 
clusion.  Analytical  treatments  of  the  interaction  phenomena 
are  practically  nonexistent,  and  those  which  exist  are  so 
superficial  as  to  be  almost  useless.  Although  much  of  the 
available  experimental  data  on  the  interaction  phenomena  are 
excellent,  in  many  cases  data  have  obtained  under  what 
amounts  to  unspecified  conditions.  The  basic  problem  would 
seem  to  be  quite  simple:  there  aren't  enough  experimental 
data.  Although  this  lack  of  information  has  been  specifically 
referenced  for  the  free-molecular- flow  regime,  it  is  apparent 
that  similar  information  is  needed  as  a  basis  for  analytical 
treatments  of  the  slip  and  transition  regimes. 

There  are  few  analytical  treatments  of  the  slip  and  trans¬ 
ition  flow  regimes.  Attempts  to  analyze  these  regimes  probably 
suffer  as  much  from  a  lack  of  a  solid  theoretical  basis  as  they 
do  from  a  lack  of  highly  accurate  data  against  which  to  test 
new  theoretical  treatments.  The  analytical  solutions  to  prob¬ 
lems  in  the  free-molecular- flow  regime  are  excellent  if  the 
appropriate  interaction  data  are  available.  Additional  work 
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is  required  on  problems  of  non-equilibrium  flows  with  non- 
Maxwellian  velocity  distributions. 

Conclusions 

At  the  present  time  there  is  a  need  for  additional  ex¬ 
perimental  data  on  the  slip,  transition  and  free-molecular- 
flow  regimes  of  rarefied  gas  dynamics.  The  need  is  especially 
great  for  conditions  corresponding  to  hypervelocity  conditions. 
In  order  to  run  controlled  experiments  in  these  regions,  it 
will  be  necessary  to  produce  certain  ranges  of  Knudsen  numbers 
and  molecular  speed  ratios,  and  to  control  the  surface  con¬ 
ditions  of  any  models  used  in  the  experimental  test  facility. 

An  accurate  knowledge  of  the  state  of  the  flow  approaching  the 
model  will  also  be  required,  and  this  may  be  one  of  the  major 
stumbling  blocks  in  obtaining  controlled  experimental  data  for 
use  as  design  information  and  as  a  basis  for  the  confirmation 
or  denial  of  new  theoretical  considerations  of  rarefied  gas 
dynamics . 
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THE  LOW  DENSITY  WIND  TUNNEl 


Configuration 

Figure  2  shows  schematically  the  low-density  wind  tunnel 
being  considered.  Gas  obtained  from  either  storage  tanks  or 
a  compressor  Is  supplied  to  a  flow  meter  and  then  passes 
through  a  heating  section,  a  settling  chamber,  and  to  the 
nozzle  section  of  the  wind  tunnel.  The  heater  will  probably 
be  of  the  arc  type.  It  is  felt  that  a  settling  chamber  of 
some  kind  is  required  for  low- density  wind  tunnels  in  order 
to  insure  that  equilibrium  has  been  achieved  prior  to  the 
nozzle  expansion  process,  and  consequently,  that  a  known  set 
of  stagnation  conditions  exist  for  the  wind  tunnel.  The  re¬ 
quirement  for  a  state  of  equilibrium  prior  to  entry  into  the 
nozzle  is  particularly  important  for  low-density  wind  tunnels, 
since  the  ranges  of  stagnation  pressures  and  stagnation  tem¬ 
peratures  which  are  encountered  inevitably  lead  to  a  combina¬ 
tion  of  equilibrium,  non-equilibrium  and  frozen  flow  in  the 
nozzle.  The  particular  combination  of  these  types  of  flow 
exerts  a  very  strong  influence  on  the  final  test  section 
conditions.  The  use  of  a  nozzle  supplied  immediately  from 
an  arc- type  of  heater  will  probably  lead  to  unknown  test- 
section  conditions,  since  the  flow  in  this  case  certainly  does 
not  achieve  an  equilibrium  condition  prior  to  entry  into  the 
nozzle,  and  therefore  the  expansion  process  in  the  nozzle  is 
unknown.  The  specific  model  to  be  tested  is  suspended  in 
the  jet  at  the  exJLt  plane  of  the  nozzle.  Upon  leaving  the 
nozzle,  the  jet  expands  into  an  additional  chamber  and  then 
in  turn  passes  through  a  precooler  section  and  finally  into 
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the  cryogenic  pumping  unit.  It  is  felt  that  a  cryogenic  pump¬ 
ing  unit  will  be  required,  since  in  order  to  achieve  flow  con¬ 
ditions  in  the  free -molecular- flow  region,  test  section 
pressures  of  the  order  of  one  micron  will  be  required.  Con¬ 
ventional  mechanical  pumping  units  are  not  suited  for  the 
volume  flow  rates  which  will  be  encountered  at  these  pressures. 
The  capacity  of  a  cryogenic  pumping  unit  is  low,  of  the  order 
of  a  few  thousand  watts,  and  the  precooler  is  a  critical  com¬ 
ponent  in  the  system  since  it  must  remove  a  sufficient  amount 
of  heat  from  the  gas  stream  so  that  the  cryogenic  unit  is  not 
overloaded.  Adequate  radiation  shielding  must  also  be  supplied 
around  the  cryogenic  unit  to  keep  the  heat  load  at  a  minimum. 
Pressure  drop  through  the  precooler  section  obviously  must  not 
be  excessive. 

Static  Pressure  Requirements 

The  requirements  for  low-density  wind  tunnels  are  identi¬ 
cal  with  the  requirements  which  have  always  been  imposed  on 
wind  tunnels;  namely,  that  a  uniform  region  of  flow  at  known 
conditions  must  be  produced.  As  has  been  mentioned  before, 
an  exact  knowledge  of  known  conditions  may  be  difficult  to 
obtain  due  to  the  instrumentation  difficulties.  Some  of  the 
quantitative  requirements  for  such  wind  tunnels  can  be  cal¬ 
culated  as  follows. 

There  exists  a  relatively  simple  relationship  between  the 
Knudsen  number,  Mach  number,  and  Reynolds  number,  vrtiich  will 
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permit  the  calculation  of  required  teat-section  pressures.  The 
relationship  Is  derived  In  Appendix  A  and  Is  presented  below: 


The  quantity  t  In  Eq.  (5)  Is  a  speed-of- sound  parameter  defined 
by  the  following  equation: 


T 


2 

a  p 
P 


(6) 


For  flow  of  a  gas  under  conditions  such  that  the  composition  Is 
frozen,  the  speed  of  sound  parameter,  t,  corresponds  to  the  is- 
entroplc  exponent  or  the  ratio  of  specific  heats,  Y  •  For  flow 
conditions  under  which  the  composition  is  controlled  by  chemical 
equilibrium  considerations,  t  and  JT  are  not  identical.  Figure 
6  of  Ref.  (5)  shows  the  speed-of-sound  parameter  for  air  under 
equilibrium  conditions  for  a  wide  range  of  pressures  and  tem¬ 
peratures.  The  compressibility  factor,  Z,  is  defined  as: 


Z 


(7) 


Values  of  Z  for  air  under  equilibrium  conditions  are  shown  in 
Fig.  1  of  Ref.  (5).  For  conditions  such  that  t  and  are  equal, 
Eq.  (5)  reduces  to  the  following  form: 


As  Y  varies  between  1.0  and  1.667,  the  coefficient  of  Eq. 

(8)  varies  between  1.277  and  1.648.  The  important  significance 
of  the  above  equations  is  that  selected  values  of  Knudsen  number 
and  Mach  number  determine  a  value  for  the  Reynolds  number  if  a 
value  of  either  t  or  Y  is  known  for  the  particular  flow  conditions. 
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It  is  interesting  to  rearrange  Eq.  (8)  into  the  following 


form: 


(9) 


where  Re  ^ 
as: 


represents  a  mean- free -path  Reynolds  number  defined 


Re 


A 


■  Kn  .  Re 


(10) 


Equation  (10)  indicates  that  for  a  desired  value  of  Mach  number 
there  is  a  single  value  of  the  mean-free-path  Reynolds  number 
which  is  required.  Figure  3  shows  Eq.  (10)  for  the  limiting 
values  of  ¥  .  In  the  range  of  Mach  numbers  from  0  to  20,  the 
mean-free-path  Reynolds  number  varies  between  0  and  a  maximum 
of  approximately  35.  Since  Eq.  (10)  shows  that  the  mean-free- 
path  Reynolds  number  is  equal  to  the  product  of  the  Knudsen 
number  and  the  usual  Reynolds  number,  it  is  easy  to  appreciate 
that  in  order  to  obtain  Knudsen  numbers  of  the  order  of  10  to 
100,  Reynolds  numbers  of  unity  or  less  will  be  encountered. 

The  Reynolds  number  is  related  to  test-section  pressure 
by  cue  following  equation: 


This  equation  is  derived  in  Appendix  A.  For  prescribed  values 
of  the  Reynolds  number  per  unit  length,  Mach  number,  and  test- 
section  static  temperature,  it  is  possible  to  calculate  the 
required  test-section  static  pressure  from  the  above  equation. 
Two  additional  assumptions  were  made  in  the  calculation  of  the 
static  pressure.  First,  the  viscosity  of  the  gas  mixture  was 
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calculated  from  the  Sutherland  equation  corresponding  to  undis¬ 
sociated  air  (or  diatomic  molecules) 


H  -  0.073225  x  10"5  - (12) 

1  +  T 

where  H  has  the  units  of  lbm/ft-sec  and  T  is  in  Rankine  tempera¬ 
ture  units.  This  approximation  has  been  shown  to  be  valid  by 
Hansen  (5) .  Second,  it  has  been  assumed  that  the  expansion  be¬ 
tween  the  stagnation  section  ahead  of  th<  nozzle  and  the  test 
section  could  be  approximated  by  an  isentropic  expansion  with  a 
constant  value  of  X  •  Values  of  Y  of  1.2,  1.4  and  1.67  were 
used  and  represent  reasonable  values  for  the  equilibrium  flow 
or  frozen-compositional  flows  which  might  occur  in  the  nozzle. 

In  Eq.  (11)  values  of  T  were  calculated  using  the  values  of  Tq, 

M  and  Y.  The  ratio  t/Z  (or  Y/Z)  was  taken  as  unity.  This  means 
that  at  low  temperatures  (or  high  M)  where  Z  and  Y  ft  1  the 
values  of  p  shown  in  Figs.  4  and  5  will  be  somevdiat  too  high. 
Further  comments  will  be  made  on  this  point  later  in  the  report. 
Minimum  and  maximum  stagnation  temperatures  of  7000  R  and 
12000R  were  used.  Figures  4  and  5  show  the  results  of  these 
calculations  in  terms  of  test-section  static  pressure  versus 
Reynolds  number  per  unit  length.  Lines  of  constant  test-section 
Mach  number  are  also  indicated,  as  are  lines  corresponding  to 
constant  values  of  the  Knudsen  number  -  characteristic  length 
product,  or  equivalently,  the  mean-free-path,  ^  .  The  values 

of  Mach  number  and  Reynolds  number  per  unit  length  correspond¬ 
ing  to  a  constant  Kn  •  L  product  were  calculated  from  Eq.  (8) . 
The  values  of  mean-free-path  were  selected  so  that,  based  on  a 
one-inch  model  characteristic  dimension,  the  three  regimes  of 
rarefied  gas  dynamics  are  covered.  For  a  stagnation  temperature 
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of  7000  R  (Fig.  4),  the  static  pressure  range  for  mean-free- 
paths  between  0.01  and  100  inches  (considering  all  Mach  num¬ 
bers)  ranges  from  approximately  1000  microns  (1  mm  Hg)  down 
to  approximately  0.01  micron.  The  main  grid  on  Fig.  4  corre¬ 
sponds  to  a  value  of  f  of  1.4.  In  addition,  two  lines  for 
a  value  of  t  of  1.2  are  shown  for  Mach  numbers  of  1  and  10. 

The  effect  on  static  pressure  of  a  change  in  the  value  of 
Jf  near  a  Mach  number  of  unity  is  extremely  small,  while  at 
a  Mach  number  of  10  there  is  an  increase  by  almost  a  factor 
of  2  in  the  required  static  pressure  as  X  is  decreased  from 
1.4  to  1.2.  These  effects  are  to  be  expected  in  view  of  the 
fact  that  from  Eqs.  (11)  »nd  (12)  the  static  pressure  de¬ 
pendency  on  temperature  is  as  follows: 

p  T  (13) 

For  Mach  numbers  near  unity,  the  ratio  of  stagnation  to 
static  temperature  can  be  written  as  follows: 

T 

T  ^  1  d*) 

Thus,  for  a  constant  value  of  the  stagnation  temperature, 
changes  in  X  should  not  affect  values  of  the  static  pressure. 
On  the  other  hand,  for  Mach  msnbers  which  are  large  compared  to 
unity,  the  ratio  of  the  stagnation  to  static  temperature  can  be 
written  as  follows: 


T 

o 

T 


M 


(15) 
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Thus,  the  static  temperature  dependency  on  #  Involves  a 
(  Y  -  1)  relationship.  For  the  values  of  1.2  and  1.4,  the 
static  temperature  ratio  for  a  given  stagnation  temperature 
should  be  in  the  ratio  of  about  2:1.  This  checks  the  varia¬ 
tion  of  static  pressure  shown  in  Fig.  4.  Figure  5  shows  the 
static  pressure  at  a  stagnation  temperature  of  12000  R.  It 
will  be  seen  that  the  values  of  static  pressure  are  about  70 
percent  larger  than  those  shown  for  7000  R.  This  effect  is 
also  expected  as  a  result  of  Eq.  (13),  since  the  static 
pressure  for  a  given  Mach  number  and  Y  should  vary  directly 
as  the  stagnation  temperature. 

In  order  to  put  these  values  of  static  pressure  in  proper 
perspective  it  can  be  mentioned  that  present-day  low-density 
wind  tunnels  using  much  lower  stagnation  temperatures,  namely, 
100-200  F,  have  been  operated  in  a  range  of  test  section 
pressures  extending  down  to  about  10  microns.  From  Fig.  4  it 
will  be  seen  that  if  a  one  inch  model  is  placed  in  a  low  den¬ 
sity  wind  tunnel  and  a  Knudsen  number  of  100  is  desired,  test 
section  pressures  between  1  micron  and  0.01  micron  will  be 
required.  Such  an  operating  point  will  put  a  severe  load  on 
existing  cryogenic  pumping  units.  Figures  4  and  5  also  in¬ 
dicate  the  necessity  for  using  extremely  small  models  in  cases 
where  large  values  of  the  Knudsen  number  are  required.  For 
example,  if  a  pumping  system  were  available  which  could  operate 
in  the  range  from  about  5  microns  to  50  microns  and  it  was 
desired  to  get  a  Knudsen  number  of  100,  Fig.  4  indicates  that 
due  to  the  pumping  restriction,  a  mean-free-path  of  about  1 
inch  would  be  possible.  The  combination  of  a  Knudsen  number 


-15- 


of  100  and  a  mean-free-path  of  1  inch  would  require  a  character¬ 
istic  length  for  the  model  of  0.01  inches.  This  is  in  fact  the 
method  by  which  large  values  of  the  Knudsen  number  have  been 
obtained  experimentally  to  date.  Curves  for  a  value  of  Y  of 
1.67  have  also  been  included  on  Figs.  4  and  5.  Again  only  the 
limiting  Mach  numbers  of  1  and  10  have  been  used. 

Two  additional  comments  should  be  made  about  Figs.  4  and  5. 
First,  values  of  the  static  pressure  corresponding  to  a  new 
stagnation  temperature  can  be  estimated  from  Eq.  (13)  using 
either  Figs.  4  or  5.  The  grid  on  these  figures  is' displaced 
vertically  to  a  new  position  such  that  the  ratio  from  Eq.  (13) 
is  produced.  Results  for  the  limiting  values  of  stagnation 
temperature  for  this  study  have  been  shown  independently  for 
illustration  purposes.  Second,  the  information  given  in  Figs. 

4  and  5  can  be  shown  on  a  single  graph.  The  static  pressure 
is  plotted  versus  Mach  number  for  a  particular  value  of  the 
Reynolds  number  per  unit  length.  At  any  new  value  of  the 
Reynolds  number  per  unit  length  the  static  pressure  can  be 
obtained  by  a  multiplication  of  the  ratios  of  the  Reynolds 
number  per  unit  length  involved.  If  the  subscript  r  is  used 
to  indicate  the  values  shown  on  this  reference  plot  and  the 
unidentified  symbols  are  used  to  indicate  new  conditions  then 
the  new  static  pressure  can  be  obtained  as  follows: 


P  - 
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Similarly  the  effect  of  a  change  In  stagnation  temperature  may 
be  calculated  from  the  following  equation: 


P  - 


(17) 


The  last  equation  is  restricted  to  high  stagnation  temperatures 
due  to  the  approximation  upon  which  Eq.  (13)  is  based.  This 
method  of  plotting  has  been  used  to  present  the  results  of  the 
approximate  boundary  layer  analysis  for  the  nozzle  in  Figs.  13 
and  14. 


Stagnation  Pressure  Requirements 

It  has  already  been  mentioned  that  non-equilibrium  flow 
effects  will  be  encountered  due  to  the  stagnation  pressures  and 
temperatures  which  are  involved.  Since  a  range  of  stagnation 
temperature  has  already  been  selected  it  is  appropriate  that 
the  possible  range  of  stagnation  pressure  also  be  estimated. 
Based  on  the  results  shown  in  Figs.  4  and  5,  it  appears  that  a 
reasonable  "average"  test-section  pressure  for  the  wind  tunnel 
might  be  about  10  microns.  This  value  was  selected  since 
static  pressures  less  than  this  value  will  require  large  wind 
tunnels  while  static  pressures  above  this  value  will  permit 
the  use  of  reasonably  sized  wind  tunnels  but  will  necessitate 
the  use  of  small  models  such  as  are  presently  used.  Fig.  6 
shows  values  of  the  nozzle  stagnation  pressure  which  will  be 
required  for  a  static  pressure  of  10  microns  and  the  Mach 
numbers  indicated  on  the  figure.  The  values  of  Y  which  were 
used  previously  have  been  used  in  Fig.  6.  For  a  Mach  number 
of  10  and  a  value  of  Y  of  1.2,  the  required  stagnation  pressure 
is  approximately  23  atmospheres.  As  Y  is  increased  to  1.4  the 
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stagnation  pressure  decreases  to  approximately  0.56  atmospheres. 
Finally,  If  a  monatomic  gas  Is  employed  the  stagnation  pressure 
Is  approximately  0.09  atmospheres.  At  this  point  It  Is  actually 
not  known  whether  or  not  a  static  pressure  of  10  microns  can  be 
tolerated  since  no  estimate  has  been  made  of  nozzle  size,  bound¬ 
ary  layer  thickness,  etc.;  however,  the  results  shown  In  Fig.  6 
do  Indicate  that  stagnation  pressures  below  approximately  100 
atmospheres  will  be  required  for  the  operation  of  the  low-den¬ 
sity  hypervelocity  wind  tunnel. 
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approximate  boundary  layer  analysis 


Introduction 

In  general,  the  nozzle  of  any  wind  tunnel  Is  designed  by 
means  of  a  two-step  procedure.  The  first  step  Involves  a  se¬ 
lection  of  the  stagnation  pressure,  the  stagnation  temperature 
and  a  desired  exit  Mach  number  for  the  nozzle.  The  specifi¬ 
cation  of  these  quantities  permits  a  study  of  the  thermodynamic 
states  and  the  isentropic  exponent  for  the  nozzle  expansion 
process.  If  the  range  of  thermodynamic  states  Involved  is 
small,  a  constant  value  for  Y  can  be  used  to  design  the  isen¬ 
tropic  contour  for  the  nozzle.  Usually  some  distribution  of 
Mach  number  along  the  centerline  of  the  nozzle  is  used  to 
generate  the  isentropic  contour.  It  is  necessary  to  keep  the 
angle  of  divergence  of  the  isentropic  contour  small  enough 
in  order  to  avoid  separation  of  the  flow  in  the  expansion 
portion  of  the  nozzle.  On  the  other  hand,  if  the  range  of 
thermodynamic  states  is  large,  the  use  of  a  constant  value 
of  Y  may  not  be  permissible,  and  it  may  be  necessary  to  use 
a  variable  value  of  Y  to  design  the  isentropic  contour.  A 
study  by  Guentert  and  Neumann  (9)  has  shown  that  significant 
errors  may  be  introduced  by  the  use  of  a  constant  value  of  y 
in  those  cases  where  large  variations  in  Y  exist. 

The  second  step  in  the  design  process  involves  a  cal¬ 
culation  of  the  boundary  layer  thickness  along  the  nozzle. 

The  quantity  of  interest  is  the  displacement  thickness  of 
the  boundary  layer,  6  .  The  boundary  layer  growth  is  cal¬ 
culated  from  an  integral  solution  of  the  boundary  layer 
equations.  The  free-stream  conditions  for  the  boundary  layer 


-19- 


calculation  are  assumed  to  be  those  from  the  isentropic  contour 
calculation.  For  the  case  of  high-density  wind  tunnels  the 
boundary  layer  thicknesses  are  usually  small  and  represent  a 
small  percentage  increase  in  the  coordinates  of  the  isentropic 
contour.  On  the  other  hand,  for  the  low-density  wind  tunnel, 
the  boundary  layer  corrections  are  large  and  may  change  the 
Isentropic  contour  coordinates  by  large  factors.  Examples  of 
such  calculations  for  low-density  wind  tunnels  are  to  be  found 
in  References  10,  11,  12.  Thus,  the  designer  of  nozzles  for 
high-density  wind  tunnels  can,  with  a  great  degree  of  certainty, 
start  a  nozzle  design  with  assurance  that  the  boundary  layer 
corrections  will  be  practically  negligible  compared  to  the  is¬ 
entropic  contour  coordinates.  However,  the  designer  of  low- 
density  wind  tunnel  nozzles  must  face  the  possibility  that  the 
boundary  layer  corrections  will  be  of  such  a  magnitude  that, 
lAien  added  to  the  isentropic  contour,  there  will  be  no  region 
of  uniform  isentropic  flow  in  vfaich  to  place  a  test  model. 

Unless  it  is  acceptable  to  utilize  such  a  nozzle  design  it  is 
necessary  to  repeat  the  entire  process  for  a  new  set  of  con¬ 
ditions  until  a  suitable  operating  point  is  found.  In  order 
to  eliminate  this  possibility  and  in  order  to  determine  suitable 
ranges  for  the  operation  of  low-density  hypervelocity  wind 
tunnels,  the  following  approximate  boundary  layer  analysis  has 
been  developed. 

Summary  of  Existing  Information 

The  analysis  and  experimental  data  for  several  nozzles  from 
the  low  density  research  groups  at  the  University  of  Toronto  and 
the  University  of  California  are  given  in  References  10,  11,  12, 
14-19.  The  analytical  techniques  tfiich  have  been  developed  for 
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boundary  layer  calculations  for  low-density  nozzles  have  been 
highly  successful  as  will  be  seen  from  an  Inspection  of  the 
experimental  results  In  these  references.  In  particular,  the 
results  of  Reference  (10)  have  been  utilized  for  boundary 
layer  calculations  In  this  report.  In  addition,  a  correlation 
based  on  the  experimentally  determined  values  of  boundary 
layer  thickness  will  be  employed  for  the  approximate  analysis 
given  here.  One  form  of  this  correlation  for  the  boundary 
layer  thickness  Is: 


6  C  _  _ C 

r'’  fwT 

In  Eq.  (18)  C  is  a  constant  which  depends  on  the  exit  Mach 
number  of  the  nozzle.  The  length  Reynolds  number  is  based 
on  the  stagnation  viscosity  rather  than  free-stream  viscosity 
at  the  exit  plane  of  the  nozzle.  Equation  (18)  has  the  same 
form  as  the  expression  for  the  laminar  boundary  layer  thick¬ 
ness  on  a  flat  plate,  namely: 


(19) 


In  Eq.  (19),  C  is  a  constant  which  depends  on  the  free-stream 
Mach  number  and  on  the  ratio  of  the  flat  plate  surface  tem¬ 
perature  to  the  free-stream  temperature .  Actually,  all  this 
latter  ratio  controls  Is  the  property  variation  in  the  laminar 
boundary  layer  for  a  specific  set  of  conditions.  If,  in 
addition,  dissociation  of  the  gas  In  the  boundary  layer  is 
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important,  then  the  free- stream  pressure  level  would  have  to 
be  introduced  as  one  of  the  controlling  variables  for  the 
constant  C.  Also,  the  Reynolds  number  in  Eq.  (19)  is  based 
on  the  free-stream  viscosity.  Manipulation  of  Eq.  (18) 
yields  the  following  equation: 


The  suggested  equivalence  between  the  two  constants  C  and  C 
then  becomes: 


C  -  C 


(21) 


It  will  be  seen  in  Eq.  (18)  that  the  axial  distance  be¬ 
tween  the  nozzle  throat  and  the  nozzle  exit  plane  has  been  used 
in  the  Reynolds  number.  Normally,  one  would  expect  that  the 
distance  along  the  wall  of  the  nozzle  would  be  the  controlling 
length  parameter  in  determining  the  boundary  layer  thickness. 
However,  the  argument  can  be  made  that  since  the  angle  of  di¬ 
vergence  in  most  wind  tunnel  nozzles  is  kept  relatively  small 
the  difference  between  the  lengths  along  the  wall  of  the  nozzle 
and  the  nozzle  axis  will  not  be  large.  Table  II  summarizes 
some  of  the  data  upon  which  Eq.  (18)  is  based.  A  check  on  the 
validity  of  Eq.  (19)  will  be  presented  later  in  this  section. 

The  analytical  boundary  layer  correction  method  which 
suggested  the  use  of  Eq.  (18)  is  presented  in  Reference  (10) . 
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This  method  uses  a  sine-arch  velocity  profile  in  order  to  solve 
the  integral  boundary  layer  equations.  The  equation  for  the 
sine-arch  velocity  profile  is: 


V 

V. 


sin 


(22) 


Based  on  Eq.  (22)  it  can  be  shown  that  the  ratio  of  the  boundary 
layer  thickness,  6  ,  to  the  displacement  thickness,  6  ,  can  be 
expressed  as  follows: 

6*  -  (1  -  B)6  (23) 


(24) 


These  results  will  now  be  used  to  develop  an  approximate  analysis 
of  the  boundary  layer  growth  in  a  low-density  wind  tunnel  nozzle. 


Summary  of  Approximate  Boundary  Layer  Analysis 

The  nomenclature  for  the  approximate  boundary  layer  analysis 
of  low-density  wind  tunnel  nozzles  is  given  in  Fig.  7.  Figure  7 a 
shows  the  divergent  section  of  the  isentropic  nozzle  contour  which 
will  be  used.  This  divergent  section  has  a  conical  shape  with  a 

if 

total  included  angle,  0.  The  nozzle  throat  diameter  is  D  and  the 
exit  diameter  of  the  isentropic  contour  is  Dgg.  The  axial  dis¬ 
tance  between  the  nozzle  throat  and  the  nozzle  exit  section  is  L^. 

Figure  7b  shows  the  final  nozzle  contour.  The  displacement  thick¬ 
et 

ness  6  has  been  added  to  the  isentropic  contour  in  order  to  ob¬ 
tain  the  final  nozzle  contour.  The  exit  diameter  of  the  final 
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contour  Is  D  .  Since  the  boundary  layer  thickness  at  the  exit 

of  the  nozzle  is  greater  than  the  displacement  thickness  at  this 

section,  not  all  of  the  exit  diameter  is  free  from  viscous 

boundary  layer  effects.  The  diameter  of  the  isentroplc  core  at 

the  exit  of  the  final  nozzle,  D1  ,  will  be  determined  by  the 

se 

values  of  6  and  D  through  the  following  relationship: 
e 

D-D  +  26*  -  D'  +  26  (25) 

e  se  se 

In  Fig.  7b  the  line  which  is  marked  "edge  of  boundary  layer"  is 
not  meant  to  indicate  the  exact  shape  of  the  boundary  layer. 

It  will  be  seen  from  the  following  analysis  that  only  the  value 
of  the  boundary  layer  thickness  at  the  exit  of  the  nozzle  is 
important. 

It  can  be  shown  that  an  expression  for  6  can  be  written  as 
follows:  (See  Appendix  B) 


The  boundary  layer  thickness  has  been  given  as  a  fraction  of  the 
exit  radius  of  the  final  contour,  rg.  This  equation  involves 
five  unknown  quantities,  namely,  B,  C,  6,  A/A*,  and  R®Dse  •  The 
angle  6  can  be  selected  at  some  reasonable  value,  while  a 
specification  of  the  Mach  number  and  V  will  determine  values  for 
B  and  A/A*.  The  quantity  C  is  determined  by  M  and  the  heat 
transfer  or  temperature  boundary  condition  at  the  nozzle  surface. 
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Finally,  a  specification  of  a  value  for  ReDge  will  permit  the 
solution  of  Eq.  (26).  There  is  a  more  direct  approach  based 
on  a  rearranged  form  of  Eq.  (26): 


Re 


Dse 


2  2 
2B  C 


b  (T7T  ‘  l)+1 
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(27) 


Although  both  equations  involve  the  same  variables,  the  use 

of  Eq.  (27)  permits  a  specification  of  a  desired  value  of 

6/rg  .  Since  the  boundary  layer  becomes  thicker  as  the 

Reynolds  number  is  decreased,  there  is  a  minimum  value  of  the 

Reynolds  number  for  which  the  boundary  layers  just  meet  and 

close  at  the  nozzle  exit  section.  For  this  condition  there 

is  no  isentropic  core  left  in  \Aiich  to  place  a  model.  The 

value  of  Re_  for  which  this  condition  exists,  or  6/r  is 
Dse  e 

unity,  can  be  written  down  directly  from  Eq.  (27) . 


<RWmln  ^|==rj  < 

It  is  also  evident  from  Eq.  (27)  that  values  of  ReDse  can 

calculated  for  any  specified  values  of  A/r  .  An  alternate 

e 

method  is  to  use  the  following  equation: 

<RW  ■  <ReD.e)mln  •  [b  («7T  *  l)  +  X]  (29) 


Hie  actual  boundary  layer  thickness  can  be  calculated  from 
the  following  expression  (See  Appendix  B)  : 
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-  1)  +  1 
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(30) 
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It  will  be  noted  that  it  is  necessary  to  specify  the  value  of 
the  Reynolds  number  per  unit  length.  However ,  this  quantity 
is  directly  related  to  the  results  shown  in  Figs.  4  and  5  and 
therefore  can  be  related  to  a  desired  operating  point  for  the 
wind  tunnel.  The  boundary  layer  thickness  corresponding  to  a 
value  of  6/r^  of  unity  follows  from  Eq.  (30): 


(*) 
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Re/L 


(31) 


Boundary  layer  thicknesses  for  values  of  6/re  different  from 
unity  can  be  calculated  from  the  following  equation  which 
results  from  combining  Eqs.  (30)  and  (31). 


6  '  <6)»  .  1  '  [b  -  «  +  l]  <32> 

re 

Once  a  value  of  6  has  been  calculated  the  displacement  thick¬ 
ness  is  immediately  known  from  Eq.  (23).  The  exit  diameter  of 
the  isentropic  contour,  the  exit  diameter  of  the  final  contour, 
and  the  diameter  of  the  isentropic  core  region  in  the  final 
contour,  can  be  calculated  from  the  following  three  equations: 


Dse- 

Re 
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Re/L 
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Some  calculated  results  based  on  the  above  analysis  will  be 
presented  later  In  this  section  of  the  report. 

Discussion  of  Assumptions  of  Approximate  Analysis 

Before  using  the  results  of  the  above  analysis  It  Is 
worthvrtiile  to  Investigate  the  validity  of  the  assumptions  for 
the  case  of  high- stagnation  temperature  flows. 

The  use  of  a  sine-arch  velocity  profile  to  solve  the  inte¬ 
gral  boundary  layer  equations  has  been  successful  in  the 
solution  of  laminar  boundary  layer  problems  at  relatively  low 
stagnation  temperatures.  Two  results  of  the  sine-arch  velocity 
profile  assumption  will  be  considered.  First,  the  correlation 
equation,  Eq.  (18),  has  been  used  to  predict  the  boundary  layer 
thickness.  Although  the  basis  for  this  equation  is  entirely 
experimental,  there  Is  an  implication  that  the  calculated  values 
of  the  boundary  layer  thickness  from  the  sine-arch  velocity  pro¬ 
file  analysis  agree  with  the  measured  values  of  the  boundary 
layer  thickness.  Second,  the  calculated  values  of  the  boundary 
layer  thickness,  6,  are  converted  to  values  of  the  displacement 
thickness,  6*,  by  the  parameter  B  which  does  depend  on  use  of 
the  sine-arch  velocity  profile  assumption. 

Figure  8  shows  values  of  the  parameter  B  as  a  function  of 
Mach  number  and  Y.  At  values  of  the  Mach  number  near  unity  the 
values  of  B  differ  only  slightly,  whereas  at  a  Mach  number  of 
10  there  is  a  factor  of  2  between  the  values  of  B  corresponding 
to  values  of  Y  of  1.1  and  1.667. 

Figure  9  shows  a  comparison  of  the  ratio  of  boundary  layer 
thickness  to  displacement  thickness  for  the  sine-arch  velocity 
profile  assumption  and  the  theoretical  values  of  Romlg  and 
Dora  (6) .  The  calculations  of  Romlg  and  Dore  were  made  for  the 
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case  of  the  laminar  boundary  layer  on  the  flat  plate  with  zero 
free-stream  pressure  gradient.  Air  was  used  as  the  working 
fluid  and  it  was  assumed  that  thermodynamic  equilibrium  exist¬ 
ed  at  all  points  in  the  boundary  layer.  The  effects  of  dis¬ 
sociation  and  property  variation  were  included.  For  a  free- 
stream  pressure  of  one  atmosphere  and  a  free -stream  temperature 
of  800  R  the  theoretical  result  is  in  good  agreement  with  the 
sine-arch  predictions  at  a  Mach  number  of  approximately  3  and 
falls  on  a  curve  for  Y  of  1.3.  As  the  Mach  number  is  increased 
the  points  gradually  move  over  to  the  curve  for  Y  of  1.2.  Hiis 
variation  would  be  expected  for  the  following  reason.  At  a 
pressure  of  1  atmosphere  and  a  temperature  of  800  R  air  is 
relatively  undissociated  and  would  have  a  value  of  Y  of  about 
1.38.  As  the  free-stream  Mach  number  is  increased,  with  free- 
stream  static  temperature  remaining  constant,  the  variation  of 
temperature  through  the  boundary  layer  would  increase  since 
the  results  shown  in  Fig.  9  were  obtained  for  a  boundary  con¬ 
dition  at  the  plate  surface  corresponding  to  an  insulated 
plate.  The  temperature  at  the  plate  surface  will  correspond 
approximately  to  the  stagnation  temperature.  Therefore,  at  a 
Mach  number  of  about  3,  the  plate  temperature  will  be  about 
2300  R.  At  this  temperature  and  a  pressure  of  1  atmosphere 
the  value  of  Y  is  about  1.3  (See  Fig.  6,  Reference  5).  The 
average  value  of  Y  in  the  boundary  layer  must  therefore  be 
somewhere  between  1.3  and  1.38.  At  this  same  free-stream 
temperature  and  a  Mach  number  of  6,  the  adiabatic-wall  tem¬ 
perature  is  about  5000  R  and  the  corresponding  value  of  Y  at 
a  pressure  of  1  atmosphere  is  about  1.15.  Thus,  the  data 
point  of  Romig  and  Dore  at  a  Mach  number  of  6  should  fall  on 
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a  curve  for  Y  of  approximately  1.2  which  it  does  as  can  be  seen 
from  Fig.  9.  The  data  points  for  free- stream  static  tempera¬ 
tures  between  1575  R  and  7000  R  are  all  above  the  sine-arch 
prediction  for  Y  -  1.1.  This  shift  would  also  be  anticipated 
based  on  the  calculations  of  Hansen  (5).  The  general  conclusion, 
based  on  the  results  shown  in  Fig.  9,  is  that  the  use  of  the 
sine-arch  velocity  profile  to  calculate  the  values  of  the  dis¬ 
placement  thickness  from  the  boundary  layer  thickness  should 
be  accurate  to  about  10%.  Additional  comments  on  this  compari¬ 
son  will  be  made  at  the  end  of  this  section. 

The  use  of  the  boundary  layer  thickness  constant  C  will 
now  be  investigated.  Figure  10  shows  values  of  C  from  three 
different  sources.  All  values  are  based  on  theoretical  ana¬ 
lyses  of  compressible  flow  over  a  flat  plate.  Two  curves  from 
the  work  of  Van  Driest  (7)  are  shown.  The  upper  curve  corres¬ 
ponds  to  an  insulated  flat  plate  and  the  lower  curve  corres¬ 
ponds  to  the  case  where  the  plate  surface  is  cooled  to  a  tem¬ 
perature  equal  to  the  free-stream  static  temperature.  Van 
Driest* s  solutions  were  obtained  for  air  and  accounted  for 
variable  fluid  properties  although  a  constant  value  of  Y  of 
1.4  was  used.  No  considerations  of  dissociation  were  involved. 
The  data  of  Romig  and  Dore  are  also  shown.  Since  Romig  and 
Dore  did  not  present  their  results  in  terms  of  the  constant  C, 
the  boundary  layer  thicknesses  given  in  Table  1  of  Reference 
(6)  were  converted  to  a  dimensionless  form  using  the  viscosity 
equation  given  in  this  report  and  values  of  the  compressibility 
estimated  from  Reference  (5) .  The  third  set  of  data  given  in 
Fig.  10  correspond  to  the  prediction  for  C  from  Howarth  (8). 

His  result  expresses  the  boundary  layer  thickness  as  a  function 
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of  Mach  number  and  v.  His  equation  is  as  follows: 
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Howarth's  results  have  been  included  since  they  give  an  expres¬ 
sion  for  the  boundary  layer  thickness  in  closed  form. 

It  will  be  seen  that  for  a  value  of  Y  of  1.4,  the  agree¬ 
ment  between  Howarth's  results  and  the  Van  Driest  solution  for 
adiabatic-wall  conditions  is  good  up  to  a  Mach  number  of 
approximately  5.  However,  at  a  Mach  number  of  10,  there  is 
about  a  40%  difference  in  the  values  of  C.  The  data  of  Romig 
and  Dore,  for  a  free-stream  static  temperature  of  800  R  are  in 
excellent  agreement  with  the  Van  Driest  solution  for  adiabatic- 
wall  conditions.  The  data  points  fall  on  the  Van  Driest  curve 
up  to  a  Mach  number  of  6  and  the  only  other  point  at  a  Mach 
number  of  about  11  has  fallen  below  the  Van  Driest  curve.  The 
agreement  at  low  Mach  numbers  is  to  be  expected  since  the 
amount  of  dissociation  in  the  boundary  layer  is  small  and  both 
studies  took  into  account  variable  fluid  properties.  The  fact 
that  the  agreement  extends  up  to  a  Mach  number  of  6  is  inter¬ 
esting  in  itself,  for  as  has  been  seen  in  Fig.  9,  the  Romig  and 
Dore  data  for  this  free-stream  static  temperature  of  800  R  de¬ 
viated  from  the  sine-arch  prediction  for  Y  =  1.4  at  a  Mach 
number  of  3  and  had  come  into  coincidence  with  the  predictions 
for  y  ■  1.2  at  a  Mach  number  of  6.  Therefore  it  appears  that 
the  ratio  of  the  boundary  layer  thickness  to  the  displacement 
thickness  is  very  sensitive  to  an  average  value  of  Y .  The 
boundary  layer  thickness  constant  C,  on  the  other  hand,  seems 
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to  be  dependent  both  on  the  actual  range  of  fluid  properties 
involved  and  on  an  average  value  of  V.  It  is  also  interest¬ 
ing  to  note  that  for  the  case  of  a  free-stream  stagnation 
temperature  of  800  R  the  first  amounts  of  dissociation  are 
not  felt  until  a  Mach  number  of  somewhat  greater  than  6  has 
been  achieved.  This  can  be  seen  from  the  plot  of  the  com¬ 
pressibility  factor  in  Fig.  1  of  Reference  (5).  A  close  in¬ 
spection  of  the  Romig  and  Dore  data  will  show  that  as  the 
free-stream  static  temperature  is  increased  the  deviations 
from  the  Van  Driest  solution  occur  at  successively  lower 
Mach  numbers.  This  observation  is  also  in  agreement  with 
the  fact  that  at  a  higher  free-stream  static  temperature  a 
lower  Mach  number  is  required  to  produce  significantly  higher 
amounts  of  dissociation  in  the  laminar  boundary  layer.  The 
curve  for  the  Van  Driest  solution  where  the  wall  temperature 
of  the  flat  plate  is  maintained  equal  to  the  free-stream 
static  temperature  has  been  included  since  a  difference  of 
about  a  factor  of  2  is  obtained  between  the  insulated-plate 
and  the  cooled-plate  cases.  If  these  theoretical  results  are 
valid,  this  means  that  a  much  thinner  boundary  layer  can  be 
maintained  in  the  nozzle  of  the  low-density  wind  tunnel  if 
nozzle  wall  cooling  is  employed.  Such  an  approach  has  been 
justified  theoretically  by  Chuan  (13)  and  in  the  operation  of 
his  low-density  wind  tunnel  at  the  University  of  Southern  Cal¬ 
ifornia  he  is  attempting  to  verify  experimentally  the  benefits 
of  nozzle  wall  cooling. 

Equation  (21)  has  been  used  to  convert  the  values  of  C 
from  the  wind  tunnel  data  to  va lue s  of  C.  In  turn,  these  values 
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of  C  have  been  compared  Co  the  Van  Driest  values  in  order  to 
verify  the  use  of  Eq.  (19)  in  the  approximate  boundary  layer 
analysis.  The  viscosity  ratio  in  Eq.  (21)  was  calculated 
from  Eq.  (12)  using  static  temperatures  obtained  at  the  Mach 
numbers  in  Table  II.  The  stagnation  temperature  for  the  data 
of  Table  II  was  taken  as  85  F  and  a  value  of  v  of  1.4  was 
used.  The  calculated  results  are  given  in  the  last  two  col¬ 
umns  of  Table  II.  The  agreement  between  the  experimental 
values  of  C  and  the  Van  Driest  values  is  excellent  up  to  a 
Mach  number  of  4  although  there  is  scatter  in  the  experimental 
values  at  this  point.  The  data  point  at  a  Mach  number  of  6 
is  about  60  percent  higher  than  the  Van  Driest  value.  However 
based  on  the  agreement  shown  in  Table  II,  Eq.  (19)  will  be 
used  in  the  approximate  boundary  layer  analysis.  Values  of  C 
will  be  taken  from  various  sources,  and  these  sources  will  be 
indicated.  Additional  comments  on  the  values  of  C  given  in 
Table  II  will  be  made  in  the  discussion  of  Fig.  26. 

Although  the  actual  processes  which  occur  in  the  nozzle 
of  the  low  density  wind  tunnel  will  be  discussed  in  more  de¬ 
tail  in  the  next  main  section  of  this  report,  a  few  comments 
are  in  order  due  to  our  discussion  of  laminar  boundary  layer 
flows  on  a  flat  plate.  The  Van  Driest  solution  represents  a 
case  in  which  a  gas  of  constant  composition  and  variable 
fluid  properties  is  used  in  the  solution  of  the  boundary 
layer  equations.  The  Romig  and  Dore  solution  represents  a 
case  in  which  a  gas  of  variable  composition  and  variable  fluid 
properties  is  used.  In  the  latter  case,  the  composition  is 
determined  by  the  requirement  that  thermodynamic  equilibrium 
exists  throughout  the  flow  field.  Intermediate  between  these 
two  extremes  and  undoubtedly  more  complex  is  the  case  where  a 
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gas  Is  used  whose  composition  is  controlled  by  finite  rates  of 
chemical  dissociation  and  recombination  and  whose  fluid  prop¬ 
erties  are  also  variable.  For  the  ranges  of  stagnation  tem¬ 
perature  and  pressure  which  will  be  encountered  in  low-density 
wind  tunnels,  the  flow  process  in  the  nozzle  will  correspond 
to  some  combination  of  equilibrium  flow,  non-equilibrium  flow, 
and  frozen-compositional  flow.  It  is  also  to  be  anticipated 
that  the  flow  will  probably  freeze  shortly  after  the  throat  of 
the  nozzle,  and  hence  most  of  the  boundary  layer  growth  in  the 
nozzle  will  occur  under  conditions  of  constant  composition  but 
varying  fluid  properties.  The  implication  of  this  situation 
is  that  perhaps  the  solutions  of  Van  Driest  or  solutions  of 
this  type  will  be  more  realistic  than  those  of  Romig  and  Dore 
where  thermodynamic  equilibrium  has  been  assumed.  Presently 
few  solutions  for  non-equilibrium  flow  conditions  over  a  flat 
plate  exist.  Based  on  this  physical  picture,  the  results  of 
Van  Driest  have  been  employed  throughout  the  next  section  on 
the  results  of  the  approximate  boundary  layer  analysis,  al¬ 
though  the  equation  of  Howarth  has  been  employed  to  show  at 
least  order  or  magnitude  changes  in  the  results  due  to 
variation  of  v.  This  decision  was  based  on  the  fact  that 
most  of  the  expansions  to  higher  Mach  numbers  are  terminated 
at  sufficiently  low  free-stream  static  temperatures  so  as  to 
be  close  to  that  temperature  used  in  the  Van  Driest  solution. 
The  principal  factor  which  might  invalidate  the  use  of  the 
Van  Driest  solution  is  that  the  actual  composition  corresponds 
to  a  frozen  flow  situation  at  a  Mach  number  of  about  1  and  a 
static  temperature  close  to  the  stagnation  temperature  used  in 
the  wind  tunnel.  If  this  composition  is  not  radically  differ¬ 
ent  from  undissociated  air,  then  the  use  of  the  Van  Driest  or 
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Howarth  solutions  should  lead  to  fairly  accurate  results. 

The  detailed  calculations  which  were  made  later  in  this 
study  indicated  that  the  Van  Driest  values  of  C  are  too  low 
at  the  higher  Mach  numbers.  The  best  values  for  C  at  the 
present  time  are  given  in  Fig.  26. 

Results  of  Approximate  Boundary  Layer  Analysis 

Fig.  11  shows  the  minimum  exit  diameter  Reynolds  number, 
ReDseJ  (for  the  isen tropic  contour)  for  which  the  boundary 
layers  just  close  at  the  exit  plane  of  the  axial  nozzle.  These 
results  are  based  on  Eq.  (28) .  The  angle  on  the  isen tropic  con¬ 
tour  has  been  selected  at  6°  (See  Fig.  7).  This  value  was 
selected  after  inspecting  isentropic  contours  for  a  variety  of 
nozzles.  The  sine-arch  velocity  profile  was  used  to  determine 
the  value  of  B.  The  Van  Driest  values  for  C  based  on  adiabatic- 
wall  conditions  were  used  for  one  of  the  curves  shown  in  Fig. 

11,  while  the  Howarth  values  for  C  were  used  in  the  other  three 
curves.  The  minimum  Reynolds  number  for  the  Van  Driest  results 
varies  between  0  and  approximately  1600  as  the  Mach  number  is 
increased  from  1  to  10.  For  the  same  value  of  Y,  namely,  1.4, 
the  Howarth  solution  shows  an  increase  of  the  Reynolds  number 
from  0  to  approximately  2900  in  the  same  Mach  number  range. 

Thus,  at  a  Mach  number  of  10  there  is  a  difference  of  almost 

2  between  the  values  of  the  Reynolds  number.  This  is  due  to 

2 

the  fact  that  the  minimum  Reynolds  number  depends  upon  C  and 
there  is  approximately  a  40%  difference  in  the  values  of  C 
given  by  these  two  solutions. 

At  a  given  Mach  number  any  value  of  Rejjse  greater  than  that 
shown  on  a  particular  curve  in  Fig.  11  will  produce  some  region 
of  isentropic  flow  at  the  exit  section  of  the  actual  nozzle. 
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Figure  12  shows  values  of  Re^se  which  are  required  to  produce 
a  given  size  isentropic  flow  area.  The  results  in  Fig.  12 
have  been  restricted  to  a  value  of  Y  of  1.4  and  a  nozzle  angle 
of  6° .  The  Van  Driest  results  for  C  have  been  used  and  the 
sine-arch  profile  approximation  for  B  has  also  been  employed. 
The  strong  variation  of  ReDge  with  6/r  at  a  given  Mach 
number  shows  that  a  rather  severe  restriction  is  placed  on  the 
wind  tunnel  designer  according  to  the  amount  of  isentropic 
core  area  which  is  desired. 

Figures  13  and  14  show  the  calculated  results  for  some 
specific  wind  tunnel  conditions  based  on  the  approximate 
boundary  layer  analysis.  The  method  of  calculating  these 
results  was  outlined  previously  in  this  section.  Each  figure 
shows  a  group  of  three  quantities  which  are  related  in  a 
rational  manner.  In  Fig.  13  the  exit  diameter  of  the  isen¬ 
tropic  contour  is  shown  first,  then  the  displacement  thick¬ 
ness  which  must  be  added  to  the  isentropic  contour  and  finally 
the  actual  exit  diameter  for  the  nozzle.  This  procedure  is 
represented  by  Eq.  (34) .  In  Fig.  14  the  three  quantities  which 
are  shown  are  the  actual  exit  diameter,  the  boundary  layer 
thickness,  and  the  diameter  of  the  isentropic  core.  These 
three  quantities  are  related  by  the  operation  represented  by 
Eq.  (35). 

The  assumptions  for  these  calculations  are  indicated  on 
each  figure.  The  only  new  condition  which  must  be  mentioned 
is  the  selection  of  a  specific  value  for  the  Reynolds  number 

3 

per  unit  length.  The  value  for  the  figures  was  taken  as  10 
inches”*.  The  form  of  presentation  in  Figs.  13  and  14  was 
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discussed  previously  relative  to  the  presentation  of  the  re¬ 
quired  static  pressures  for  specific  test  section  conditions, 
namely  Figs.  4  and  5.  This  new  form  has  been  used  since  the 
number  of  curves  which  were  to  be  presented  would  have  pre¬ 
sented  an  unintelligible  cross-hatching  if  done  in  the  manner 
similar  to  Figs.  4  and  5.  The  numerical  values  given  Figs. 

13  and  14  can  be  changed  as  a  function  of  the  value  of  the 
Reynolds  number  per  unit  length  by  simple  proportional  scal¬ 
ing.  If  the  Reynolds  number  per  unit  length  is  increased  by 
a  factor  of  10,  all  diameters  are  reduced  by  a  factor  of  10, 
etc.  The  fact  that  a  proportional  scaling  with  Reynolds 

number  per  unit  length  is  involved  indicates  why  the  partic- 
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ular  value  of  10  inches  was  used.  This  can  be  seen  from 
the  values  of  the  actual  exit  diameters  which  are  encountered 
for  these  conditions.  If  Reynolds  numbers  per  unit  length  of 
either  10  or  100  inches  ^  had  been  used  the  actual  exit  dia¬ 
meters  would  have  become  very  large. 

The  exit  diameters  for  the  isentropic  contours  shown  in 

Fig.  13  vary  between  approximately  0.1  inch  and  1.4  inches 

for  the  case  where  the  boundary  layers  just  close  at  the  exit 

plane  of  the  actual  nozzle  or  a  value  of  6/r^  of  unity.  As 

6/r  is  reduced  from  unity  to  0.8,  at  Mach  number  of  10,  the 
6 

exit  diameter  for  the  isentropic  contour  increases  from  about 
1.4  to  approximately  7  inches.  This  factor  of  5  was  encount¬ 
ered  earlier  in  Fig.  12  where  the  minimum  Reynolds  number  for 
the  isentropic  contour  was  presented.  The  displacement  thick¬ 
nesses  which  must  be  added  to  the  isentropic  contour  dimen¬ 
sions  are  shown  in  the  middle  section  of  Fig.  13.  At  a  Mach 
number  of  10  and  a  value  of  6/rg  of  unity  the  displacement 
thickness  is  about  2.9  inches.  This  means  that  a  radial  dis- 
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placement  thickness  of  2.9  inches  must  be  added  to  the  exit 

radius  of  the  isentropic  contour  which  is  approximately  0.7 

inches.  Thus,  the  displacement  thickness  is  approximately 

four  times  the  exit  radius  of  the  isentropic  contour.  For 

a  Mach  number  of  10  and  a  value  of  6/r  of  0.8  the  exit  radius 

e 

of  the  isentropic  contour  is  about  3%  inches  while  the  dis¬ 
placement  thickness  is  approximately  6%  inches.  The  actual 
exit  diameter  is  about  20  inches.  Thus,  for  these  latter 
conditions,  the  displacement  thickness  is  about  twice  the 
exit  radius  of  the  isentropic  contour.  These  figures  again 
illustrate  that  the  boundary  layers  which  are  encountered  in 
low-density  wind  tunnel  nozzles  are  appreciable  compared  with 
either  the  isentropic  contour  dimensions  or  the  actual  nozzle 
contour  dimensions.  These  actual  exit  diameters  are  reason¬ 
able  and  represent  nozzle  dimensions  which  could  ve  used 
without  overloading  either  the  heater  or  pumping  capacities 
of  most  hypersonic  flow  installations. 

Figure  14  shows  the  actual  exit  diameter  of  the  nozzle, 
the  boundary  layer  thickness  and  the  isentropic  core  diameter 
at  the  exit  plane  of  the  actual  nozzle.  The  actual  exit  dia¬ 
meters  are  those  shown  in  Fig.  13.  The  values  of  the  boundary 
layer  thickness  follow  immediately  from  the  actual  exit  dia¬ 
meters  since  the  curves  shown  in  the  middle  section  of  Fig.  14 
correspond  to  specific  values  of  6/r^.  However,  in  the  actual 
calculation  of  the  results  shown  in  Fig.  14  the  value  of  6  was 
calculated  from  Eq.  (30)  and  a  check  was  made  to  see  that  the 
assumed  value  of  6/rg  was  satisfied.  The  isentropic  core  dia¬ 
meter  is  shown  in  the  third  section  of  Fig.  14  and  the  values 
for  6/re  of  unity  are  obviously  zero.  For  a  value  of  6/r^  of 
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0.9  the  isentropic  core  diameter  varys  from  approximately  0.03 
inch  to  about  1.2  inches  as  the  Mach  number  is  varied  between 
2  and  10.  For  a  higher  value  of  6/rg,  namely,  0.8,  the  isen¬ 
tropic  core  diameter  varys  between  0.09  inch  and  about  4  inches 
over  the  same  Mach  number  range.  The  size  of  the  isentropic 
core  is  certainly  large  enough  for  models  of  reasonable  size. 

Conclusions 

The  results  of  the  approximate  boundary  layer  analysis 
shown  in  Figs.  13  and  14  can  now  be  used  to  draw  some  conclu¬ 
sions  about  possible  operating  ranges  and  test  section  conditions 
for  low  density  wind  tunnels.  It  is  obvious  that  no  where  in 
the  boundary  layer  analysis  has  specific  mention  been  made  of 
nozzle  stagnation  pressures  or  temperatures.  These  numbers 
actually  enter  the  analysis  only  as  they  effect  the  values  of 
V  and  C  which  are  used  in  the  boundary  layer  analysis  and 
thirdly,  as  they  effect  the  value  of  Reynolds  number  per  unit 
length.  It  will  be  seen  in  the  next  main  section  of  this  re¬ 
port  that  due  to  the  combination  of  equilibrium,  non-equilibrium 
and  frozen  compositional  equilibrium  flows  which  are  accounted 
in  nozzles  of  low  density  wind  tunnels  that  a  value  of  v  of  1.3 
perhaps  should  have  been  employed.  However,  conclusions  based 
on  a  value  of  Y  of  1.4  will  not  be  significantly  different  from 
those  for  a  value  of  v  of  1.3.  Since  the  results  shown  in  Figs. 

13  and  14  have  been  based  on  a  value  of  the  Reynolds  number  per 
3  -1 

unit  length  of  10  inches  it  is  possible  to  determine  what  the 
allowable  operating  range  for  the  low  density  wind  tunnel  will 
be.  In  other  words,  the  Reynolds  number  per  unit  length  is  the 
important  link  between  the  approximate  boundary  layer  analysis 
and  the  operating  conditions  as  given  in  Figs.  4  and  5  or  charts 
similar  to  these. 
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In  Fig.  4  for  a  stagnation  temperature  of  7000  R  and 
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Reynolds  number  per  unit  length  of  10  inches  ,  it  will  be 
seen  that  a  Mach  number  of  10  the  mean-free-path  is  about 
0.02  inch  and  the  static  pressure  at  the  nozzle  exit  is  about 
80  microns.  Thus  a  model  having  a  characteristic  length  of 
0.02  inch  would  have  a  Knuds en  number  of  1  under  these  con¬ 
ditions.  If  the  characteristic  length  was  reduced  to  0.002 
inch,  the  Knudsen  number  would  increase  to  10.  Both  of  these 
Knudsen  numbers  are  in  the  transitional  flow  regime  but  the 
model  size  cannot  be  called  large.  In  fact,  these  model  sizes 
are  comparable  to  currently  used  sizes  in  low-density  wind 
tunnels . 

From  Fig.  14  it  will  be  seen  that  if  the  exit  diameter 

of  the  nozzle  is  about  12  inches,  corresponding  to  the  value 

of  6/rg  of  0.9,  then  the  isen tropic  core  region  will  have  a 

diameter  of  approximately  1.3  inches.  Thus,  both  of  the  model 

sizes  mentioned  could  be  placed  in  a  region  of  uniform  flow  at 

isentropic  flow  conditions.  Obviously,  if  an  exit  diameter  for 

the  actual  nozzle  of  120  inches  could  be  tolerated,  then  the 

Reynolds  number  per  unit  length  could  be  decreased  to  about  a 
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value  of  10  inches  .  This  would  increase  the  mean-free-path 
to  0.2  inch  and  the  model  size  would  be  increased  by  a  factor 
of  10  keeping  the  same  value  of  Knudsen  number  or  the  model 
size  could  be  maintained  thus  increasing  the  Knudsen  numbers 
by  the  same  factor  of  10.  The  isentropic  core  region  would 
have  a  diameter  of  some  13  inches.  Since  the  size  of  the  is¬ 
entropic  core  region  is  somewhat  large  compared  to  the  allow¬ 
able  model  sizes  it  is  then  possible  to  consider  a  reduction  in 
the  value  of  &/r^  that  is  required.  The  practical  implication 
of  such  a  reduction  is  that  the  exit  diameter  of  some  120 
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inches  can  probably  be  brought  down  to  a  more  realistic  value 

without  changing  either  the  Reynolds  number  per  unit  length 
2  -1 

from  10  inches  or  increasing  the  mean-free-path  from  0.2 
inch. 

It  is  also  possible  to  speculate  on  the  effects  of  un¬ 
certainties  in  the  approximate  boundary  layer  analysis  or  in 
certain  methods  of  boundary  layer  reduction.  For  example,  a 
curve  was  shown  in  Fig.  10  for  C  corresponding  to  a  cooled 
nozzle  wall.  This  curve  was  lower  by  a  factor  of  approximately 
2  than  the  curve  for  C  corresponding  to  the  Van  Driest  adia¬ 
batic-wall  conditions.  Since  C  centers  the  expressions  for  the 
minimum  value  of  Re^ge  and  the  boundary  layer  thickness  as  a 
squared  term,  the  values  shown  for  lengths  in  Figs.  13  and  14 
would  be  reduced  by  approximately  4.  It  can  be  estimated  that 
by  proper  juggling  of  model  size  and  wall  cooling  conditions 
that  the  Reynolds  number  per  unit  length  could  be  driven  down 
to  a  value  of  approximately  10  inches  ^ .  Under  such  conditions, 
at  a  stagnation  tenqperature  of  7000  R,  a  mean-free-path  of  1 
inch  would  be  possible  and  the  static  pressure  in  the  test 
section  would  be  about  one  micron.  Similarly,  for  a  stagnation 
temperature  of  12000R  (Fig.  5)  the  same  estimates  will  be  valid. 
It  does  not  appear  that  it  will  be  possible  to  test  a  model  hav¬ 
ing  a  characteristic  dimension  of  0.1  inch  in  a  low-density 
wind  tunnel  at  high  Mach  numbers  under  conditions  such  that  the 
Knudsen  number  based  on  this  length  is  100! 
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real-gas  effects  in  low-density  wind  tunnels 

Description  of  Real-Gas  Effects 

There  are  two  limiting  processes  by  which  the  character¬ 
istics  of  a  nozzle  expansion  process  can  be  calculated.  For 
states  corresponding  to  a  high  pressure  and  low  temperature, 
the  gas  composition  is  constant  during  the  expansion  process 
and  it  is  only  necessary  to  account  for  variations  of  various 
thermodynamic  properties  with  temperature.  This  approximation 
to  a  real  expansion  process  is  called  the  frozen-flow  approxi¬ 
mation.  It  should  be  noted  that  the  reference  to  a  frozen 
condition  refers  specifically  to  the  composition  of  a  gas. 

This  condition  may  be  accompanied  by  either  variable  or  con¬ 
stant  thermodynamic  properties.  For  states  corresponding  to 
a  low  pressure  and  a  high  temperature  the  gas  composition  may 
no  longer  remain  constant  throughout  the  expansion  process. 

If  the  expansion  process  is  a  "gradual"  one,  then  at  each 
section  in  the  nozzle  a  condition  of  thermodynamic  equilib¬ 
rium  will  exist.  This  approximation  to  a  real  expansion  pro¬ 
cess  is  called  the  equilibrium- flow  approximation.  In 
addition  to  variable  gas  composition  throughout  the  nozzle, 
the  thermodynamic  properties  will  also  be  variable.  Both  of 
these  limiting  approximations  are  isentropic,  and  the  nozzle 
expansion  process  may  be  calculated  independently  of  a 
specific  nozzle  contour.  For  the  ranges  of  stagnation  pressure 
and  stagnation  temperature  encountered  in  low-density  wind 
tunnels,  it  is  found  that  neither  of  the  limiting  approxi¬ 
mations  gives  an  accurate  picture  of  the  nozzle  expansion 
process.  Consideration  must  be  given  to  the  actual  rates  at 
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which  the  molecular  and  atomic  species  of  the  gas  are  dissociat¬ 
ing  and  recombining. 

Figure  15  shows  schematically  the  composition  and  static 
temperature  variations  in  a  hypersonic  nozzle.  It  is  assumed 
that  a  pure  diatomic  gas  is  used.  The  quantity  a  is  a  local 
mass  fraction  of  gas  dissociated  and  defined  as 


a 


_ (AJ  . 

(A)  +  2(A2) 


(37) 


where  A  is  the  atomic  concentration  and  is  the  molecular  con¬ 
centration.  Both  a  and  the  static  temperature  ratio  are  plotted 
as  a  function  of  area  ratio  based  on  the  throat  area  of  the 
nozzle.  Negative  values  of  the  area  ratio  are  upstream  of  the 
throat  section  and  positive  areas  are  in  the  divergent  section 
of  the  nozzle.  The  local  mass  fraction  of  gas  dissociated  at 
stagnation  conditions  is  a^.  For  an  expansion  process  in  which 
the  chemical  dissociation  and  recombination  rates  are  zero  a 
remains  constant  along  the  length  of  the  nozzle  at  the  stagnation 
state  value.  This  is  shown  by  the  upper  curve  in  Fig.  15a.  If 
the  nozzle  expansion  process  occurs  under  thermodynamic  equili¬ 
brium  conditions,  then  a.  follows  the  lower  curve  in  Fig.  15a 
labeled  "infinite  reaction  rate".  Since  under  the  particular 
conditions  involved  for  this  study  neither  of  these  approximations 
will  be  valid,  one  must  consider  the  third  curve  shown  in  Fig.  15a 
As  the  gas  passes  through  the  convergent  section  of  the  nozzle, 
the  expansion  rate  is  low  and  the  curve  for  the  actual  expansion 
follows  closely  the  curve  for  the  equilibrium- flow  approximation. 
At  some  point,  A  on  the  equilibrium  curve  the  expansion  process 
becomes  sufficiently  rapid  so  that  the  actual  rates  of  dissoc¬ 
iation  and  recombination  become  controlling.  The  balance  of 
these  rate  processes  is  such  as  to  produce  a  value  of  a  at  any 
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section  of  the  nozzle  which  is  intermediate  between  the  frozen- 
flow  and  the  equilibrium  flow  values.  As  the  nozzle  expansion 
process  proceeds,  the  Point  B  is  reached  such  that  the  com¬ 
position  of  the  gas  becomes  frozen.  The  expansion  process  be¬ 
yond  Point  B  may  be  calculated  by  a  frozen-flow  approximation. 
If  the  actual  nozzle  exit  conditions  for  the  nozzle  are  located 
downstream  of  Section  B,  it  is  necessary  to  consider  three  dis¬ 
tinct  processes  in  the  nozzle.  The  process  from  the  stagnation 
state  to  Point  A  is  an  infinite-reaction  rate  or  equilibrium- 
flow  process,  the  process  from  Point  A  to  Point  B  is  a  finite- 
reaction  rate  or  non-equilibrium- flow  process,  aid  the  region 
from  Point  B  to  the  nozzle  exit  plane  is  a  zero-reaction  rate 
or  frozen-compositional  flow  process. 

Since  the  static  temperature  is  the  flow  parameter  most 
affected  by  the  non-equilibrium- flow  conditions,  it  is  shown 
schematically  in  the  lower  half  of  Fig.  15.  The  upper  curve 
in  Fig.  15b  corresponds  to  the  equilibrium- flow  process.  The 
lower  curve  in  Fig.  15b  corresponds  to  the  frozen-flow  process 
starting  at  stagnation  conditions.  Two  intermediate  curves 
have  been  shown  in  Fig.  15b.  The  lower  curve  of  this  pair 
corresponds  to  the  actual  non-equilibrium- flow  process  and  it 
deviates  from  the  equilibrium  curve  at  Point  A  as  did  the 
curve  for  a  in  the  upper  half  of  Fig.  15.  The  remaining  curve 
in  Fig.  15b  represents  an  approximation  which  is  often  used  to 
the  actual  non-equilibrium  expansion  process.  A  freezing  cri¬ 
terion  is  selected  such  that  based  on  the  equilibrium  expansion 
curve  which  may  be  easily  obtained,  it  is  possible  to  predict 
when  "appreciable"  freezing  has  occurred.  This  defines  a 
Point  C  on  the  equilibrium  curve  and  some  value  of  the  fraction 
dissociated.  The  expansion  process  downstream  of  Point  C  is 
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assumed  to  be  s  frozen  compositional  flow  (corresponding  to  a  ) 

c 

and  account  Is  usually  taken  of  variable  thermodynamic  prop¬ 
erties  during  the  remainder  of  the  expansion  process  from  Point 
C.  the  fourth  curve  In  Fig.  15b  corresponds  to  this  approxi¬ 
mation  to  the  actual  non-equillbrlum- flow  process.  The  close¬ 
ness  of  the  curve  for  the  approximate  process  to  the  curve  for 
the  actual  non- equilibrium- flow  process  will  obviously  depend 
on  the  accuracy  of  the  freezing  criterion,  on  the  size  of  the 

Increment  between  a  and  a  ,  and  on  the  area  ratios  between 

a  B 

these  two  points.  The  theoretical  treatments,  calculated 
results,  and  experimental  data  on  the  non-equilibrium- flow  pro¬ 
cesses  have  been  presented  in  References  20  through  48. 

Typical  Results 

Figures  16,  17,  18  and  19  show  the  results  of  calculations 
for  non-equilibrium  expansion  processes  of  oxygen  In  a  15° 
hyperbolic  nozzle  (  Reference  35  ), 

The  first  three  figures  show  distribution  of  composition  along 
the  nozzle.  The  stagnation  temperature  for  these  calculations 
was  5900  K,  or  approximately  10600  R.  Fig.  16  shows  the  results 
for  a  stagnation  pressure  of  approximately  1  atm.  For  these 
conditions  the  oxygen  is  almost  completely  dissociated  and  aQ  is 
0.9937.  The  difference  between  the  results  for  the  finite- 
reaction-rate  calculations  and  the  infinite-reaction-rate  cal¬ 
culations  are  quite  large.  The  flow  freezes  very  shortly  down¬ 
stream  of  the  throat  at  an  area  ratio  of  approximately  1.4. 

After  freezing  has  occurred,  the  value  of  a  is  approximately 
0.985.  On  the  other  hand,  for  the  equilibrium  expansion  the 
value  of  a  has  dropped  to  0.86  at  an  area  ratio  of  approximately 
2.  In  Fig.  17  where  the  stagnation  pressure  has  been  raised  to 
9.4  atm.  about  94  percent  of  the  oxygen  has  been  dissociated  at 
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stagnation  conditions.  For  these  conditions  freezing  occurs  at 
an  area  ratio  of  approximately  3,  and  the  value  of  a  at  freezing 
is  approximately  0.85.  Finally  in  Fig.  18  where  the  upstream 
stagnation  pressure  has  been  increased  to  82  atm.  the  fraction 
dissociated  at  stagnation  conditions  is  about  69  percent.  Freez¬ 
ing  does  not  occur  until  an  area  ratio  of  about  10  is  reached 
and  the  frozen  composition  corresponds  to  a  value  of  a  of  about 
0.45.  Thus  it  can  be  seen  that  the  position  of  freezing  and  the 
composition  at  the  freezing  point  are  strong  function  of  stag¬ 
nation  pressure  at  this  temperature  of  5900  K. 

The  static  temperature  distributions  shown  in  Fig.  19  are 
similar  to  those  presented  schematically  in  Fig.  15b.  However, 
the  difference  in  the  final  static  temperature  at  a  given  area 
ratio  is  illustrated  dramatically  in  a  quantitative  manner  in 
Fig.  19.  At  an  area  ratio  of  1000,  for  example,  the  static  tem¬ 
perature  ratio  for  the  equilibrium  or  infinite-reaction-rate  ex¬ 
pansion  is  about  0.45  while  that  for  the  frozen  expansion  pro¬ 
cess  is  only  0.05.  If  the  nozzle  expansion  process  is  taken  to 
a  sufficiently  large  area  ratio,  it  will  be  found  that  the 
difference  in  the  velocity  for  the  actual  nozzle  expansion  pro¬ 
cess  and  the  equilibrium  expansion  process  is  not  large.  How¬ 
ever,  the  strong  dependency  of  static  temperature  on  the  partic¬ 
ular  expansion  process  will  produce  a  greatly  different  Mach 
number  at  a  given  area  ratio  due  to  the  dependency  of  the  speed 
of  sound  on  temperature.  The  pressures  and  densities  are  less 
dependent  on  the  particular  expansion  process  than  is  the  tem¬ 
perature. 

The  preceding  results  for  the  non-equilibrium  flow  of  ox¬ 
ygen  in  a  hypersonic  nozzle  are  interesting,  since  they  show  in 
a  quantitative  manner  the  nature  of  the  non-equilibrium  flow 
process.  However,  results  are  available  for  the  non-equilibrium 
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flow  of  air  in  a  hypersonic  nozzle  (20).  Some  of  these  results 
are  presented  below.  The  results  for  a  stagnation  temperature 
of  6000  K  and  a  stagnation  pressure  of  100  atm.  have  been  se¬ 
lected  due  to  their  closeness  to  the  desired  stagnation  conditions 
for  low-density  wind  tunnels.  The  nozzle  shape  Is  hyperbolic. 
Although  no  curves  for  the  composition  of  the  flow  are  presented, 
these  results  are  given  In  a  tabular  and  graphical  manner  In 
Reference  (20) .  Figure  20  shows  the  static  temperature  ratio 
versus  area  ratio  for  these  conditions.  Curves  for  the  equili¬ 
brium-flow  and  frozen- flow  approximations  are  shown.  For  the 
frozen- flow  curve  the  flow  is  assumed  to  be  frozen  at  the 
throat  composition  conditions.  Two  intermediate  curves  are 
shown,  corresponding  to  two  calculated  points  for  freezing  in 
the  nozzle.  The  relaxation  length  criteria  of  Reference  (35) 
has  been  used  to  determine  several  freezing  points  along  the 
equilibrium-expansion  curve.  The  quantity  1  is  a  length  which 
characterizes  the  nozzle  geometry  and  is  the  ratio  of  the  half- 
height  of  the  nozzle  at  the  throat  section  (in  centimeters) 
divided  by  the  tangent  of  the  nozzle  semi-angle  (See  Reference 
(20)).  The  freezing  point  is  shown  for  two  assumed  values  for 
1.  The  curves  in  Fig.  20  are  similar  to  those  shown  in  Fig. 

15b  and  Fig.  19.  The  point  of  departure  for  either  frozen  curve 
from  the  equilibrium  curve  corresponds  to  Point  C  in  Fig.  15a  or 
15b.  Once  more  it  will  be  seen  that  at  a  particular  value  of  the 
area  ratio  the  value  of  the  static  temperature  is  greatly  affect¬ 
ed  by  the  prior  history  of  the  flow  before  it  reaches  a  given 
section.  Since  most  of  our  past  experience  with  nozzles  expan¬ 
sion  processes  has  been  based  on  the  concept  of  some  average 
value  of  the  isentropic  exponent  for  the  overall  process,  three 
lines  of  constant  Y  have  been  imposed  on  Fig.  20.  Curves  for 
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values  of  Y  of  1.2,  1.3  and  1.4  are  shown.  For  nozzle  sections 
near  the  throat  It  will  be  seen  that  the  equilibrium  expansion 
curves  follow  a  curve  for  a  value  of  Y  of  about  1.2.  As  the 
expansion  proceeds  along  the  nozzle ,  the  equilibrium  curve 
first  deviates  toward  values  of  Y  less  than  1.2,  then  reverses 
itself  toward  the  curve  for  1.2,  passes  underneath  this  curve, 
and  eventually  reaches  values  of  temperature  corresponding  to 
a  value  of  Y  of  1.4.  This  behavior  of  the  curve  for  the  equi¬ 
librium  process  is  to  be  expected  in  view  of  the  values  of  Y 
given  in  References  (5)  and  (49) .  The  curve  for  the  frozen 
flow  (composition  frozen  at  throat)  falls  between  the  two  curves 
for  Y  of  1.3  and  1.4.  Near  the  throat  section  the  frozen  curve 
is  closer  to  the  1.3  curve,  tdiile  during  the  expansion  process 
throughout  the  remainder  of  the  nozzle,  the  frozen-equilibrium 
curve  gradually  shifts  over  to  the  curve  for  a  value  of  Y  of 
1.4. 

Figure  21  shows  some  additional  details  from  Fig.  20.  Two 
curves  corresponding  to  values  of  y  of  1.3  and  1.4  have  been 
drawn  starting  at  the  point  of  departure  of  the  frozen- flow 
curves  from  the  equilibrium- flow  curve.  For  both  frozen  flow 
curves  it  will  be  seen  that  near  the  point  of  departure  the 
initial  phase  of  the  frozen-flow  expansion  process  is  character¬ 
ized  by  a  value  of  Y  somewhat  less  than  1.3,  and  at  later  stages 
in  the  expansion  process  there  is  a  gradual  shift  to  a  curve 
having  a  value  of  Y  equal  to  1.4.  This  means  that  if  one 
follows  the  nozzle  expansion  process  from  the  stagnation  region 
to  the  point  of  freezing,  a  variation  in  Y  is  observed  which 
has  the  following  general  characteristics.  For  that;  phase  of 
the  expansion  process  corresponding  to  an  equilibrium  expansion. 
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the  values  of  y  are  liable  to  be  about  1.2  or  leas.  As  one  passes 
from  the  equilibrium  curve,  through  the  freezing  point  and  along 
the  frozen-flow  curve,  a  discontinuity  in  the  value  of  Y  occurs, 
such  that  there  is  a  jump  from  a  fairly  low  value  of  1.2  or  less 
to  a  value  which  is  closer  to  1.3.  Following  along  the  frozen 
flow  curve,  Y  undergoes  a  gradual  transition  from  1.3  to  1.4. 
Obviously,  these  remarks  depend  on  the  particular  set  of  nozzle 
stagnation  conditions  and  to  a  certain  extent  on  nozzle  geometry. 
The  values  of  Y  for  the  equilibrium  expansion  as  has  been  mention¬ 
ed,  are  to  be  expected  based  on  many  calculations  for  air  in 
thermodynamic  equilibrium.  The  values  of  Y  for  the  frozen  flow 
process  are  also  to  be  expected  in  view  of  the  following  ex¬ 
pressions  for  Y  given  in  Reference  (48) .  For  the  case  of  frozen 
composition  and  frozen  molecular  vibrations,  Y  should  have  a 
value  of: 
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(38) 


where  y  is  the  number  of  molecules  per  initial  molecule  or  air 
m 

and  ya  is  the  number  of  atoms  per  Initial  molecule  of  air. 
Equation  (38)  predicts  values  of  Y  between  1.4  and  1.67.  On  the 
other  hand,  if  the  composition  is  frozen  but  the  molecular  vi¬ 
brations  are  in  equilibrium,  a  value  of  Y  given  by  the  following 
equation  is  to  be  expected: 
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Equation  (39)  predicts  values  of  Y  ranging  from  approximately 
1.29  to  1.67.  Since  the  results  shown  in  Figs.  20  and  21  were 
calculated  assuming  molecular  vibrations  were  in  equilibrium, 
the  results  of  Eq.  (39)  should  apply,  and  the  observed  values 
of  y  check  those  predicted  by  Eq.  (39). 

Figure  22  shows  the  velocity  as  a  function  of  area  ratio 
for  the  same  nozzle  stagnation  conditions.  Up  to  an  area  ratio 
of  5,  at  vrtiich  point  the  velocity  is  approximately  10,500  ft/sec, 
there  is  no  difference  between  any  of  the  nozzle  flow  approxi¬ 
mation  curves.  At  area  ratios  above  this  value,  the  curves 
begin  to  diverge,  but  it  appears  that  the  maximum  difference 
between  any  of  the  curves  is  only  of  the  order  of  about  10  or 
15  percent.  Figure  23  shows  the  Mach  number  as  a  function  of 
area  ratio.  The  velocity  of  sound  for  these  curves  has  been 
calculated  from  the  usual  equation,  namely: 


The  equation  was  solved  using  either  the  frozen-flow  or  the 
equilibrium- flow  derivative.  The  important  implication  of 
Fig.  23  is  that  for  a  desired  value  of  Mach  number  a  signifi¬ 
cantly  different  area  ratio  is  required,  once  more  depending 
on  whether  the  flow  is  in  equilibrium  or  frozen.  Based  on  the 
preceding  requirements  for  low-density  wind  tunnels,  and  the 
approximate  boundary  layer  analysis,  it  is  obvious  that  the 
use  of  a  shorter  nozzle  will  produce  thinner  boundary  layers 
and  therefore  a  wider  range  of  satisfactory  operation  for  the 
wind  tunnel.  Figure  24  shows  the  static  pressure  ratio  for 
the  hypersonic  nozzle.  Factors  of  between  2  and  3  are  encoun¬ 
tered  in  the  static  pressure,  depending  on  the  actual  flow 
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process  in  the  nozzle.  Finally,  Fig.  25  shows  the  density 
ratio  for  the  nozzle,  and  it  is  found  that  the  density  is  the 
flow  parameter  which  is  least  affected  by  the  actual  nozzle 
expansion  process. 

Conclusions 

In  view  of  the  ranges  of  stagnation  pressures  and  tem¬ 
peratures  which  are  required  for  low-density  wind  tunnels,  it 
will  be  necessary  to  consider  chemical  kinetic  effects  and  the 
non-equilibrium  flow  processes  vAiich  are  produced.  Current 
calculations  of  non-equilibrium  air  flows  have  not  been  made 
at  stagnation  pressures  below  100  atm.  Based  on  the  results 
presented  earlier  in  this  report,  it  now  appears  that  these 
calculations  should  be  extended  into  the  range  of  stagnation 
pressure  between  0.1  and  100  atm.  The  range  of  stagnation 
temperature  which  should  be  considered  extends  from  3000  K  to 
10000  K.  The  detailed  curves  for  the  equilibrium,  non-equili¬ 
brium  and  frozen-flow  processes  show  that  a  wide  range  of 
values  of  Y  are  encountered;  namely,  from  about  1.1  to  1.4. 

The  upper  limit  of  1.4  was  selected  baced  on  the  calculations 
presented  in  this  report.  However,  at  higher  stagnation  tem¬ 
peratures,  where  larger  amounts  of  dissociation  and  ionization 
to  monatomic  species  will  be  encountered,  the  upper  limit  of 
Y  will  rise  towards  1.67.  This  range  of  Y  creates  somewhat  of 
a  dilemma  for  the  nozzle  designer.  In  order  to  calculate 
nozzle  exit  conditions,  a  value  of  Y  of  perhaps  1.4  may  be 
appropriate  when  used  with  the  upstream  stagnation  conditions. 
However,  to  the  nozzle  contour  designer  a  value  of  y  of  1.3 
may  be  more  characteristic  for  an  isentropic  contour  upon  which 


-51- 


to  base  the  actual  nozzle  design.  The  real  solution  to  this 
problem  is  to  design  nozzles  such  that  the  actual  thermodynamic 
processes  are  connected  directly  to  the  isentropic  contour  de¬ 
sign  and  to  any  boundary  layer  correction  method  which  will  be 
employed. 
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DETAILED  NOZZLE  CALCULATIONS 


Introduction 

Although  the  calculated  results  from  the  approximate  bound¬ 
ary  layer  analysis  permitted  a  definition  of  some  allowable 
operating  conditions  for  low-density  wind  tunnels,  detailed 
nozzle  calculations  are  required  to  check  on  the  accuracy  of 
the  approximate  analysis  and  to  provide  final  nozzle  contours 
and  system  performance  characteristics  for  the  actual  low-density 
wind  tunnels.  Three  nozzle  contours  were  designed,  including 
appropriate  boundary- layer  corrections.  The  isentropic  contours 
for  these  nozzles  were  obtained  from  two  different  sources.  Refs. 
(50)  and  (51)  .  Stagnation  temperatures  of  5500  R,  6730  R,  and 
12500  R  were  used.  The  exit  Mach  number  for  all  three  nozzles 
was  about  10.  The  boundary  layer  correction  method  of  Ref.  (10) 
was  used,  since  nozzles  designed  by  this  method  had  shown  ex¬ 
cellent  results  in  previous  experimental  tests. 

Nozzle  Design  Procedure 

Isentropic  contours  were  obtained  from  two  sources.  Ref¬ 
erence  (50)  presents  an  isentropic  contour  with  an  exit  Mach 
number  of  10.068.  In  addition  to  this  contour,  nozzle  coordinates 
were  obtained  from  the  Sandburg-Serrell  Corporation,  Ref.  (51). 

The  isentropic  contour  corresponding  to  an  exit  Mach  number  of  10 
and  a  value  of  Y  of  1.2  was  used.  The  details  of  the  boundary 
layer  correction  method  are  given  in  Ref.  (10).  Only  the  import¬ 
ant  details  will  be  sunmarized  here.  The  boundary  layer  thick¬ 
ness,  &,  corresponding  to  flow  conditions  at  any  section  along 
the  isentropic  contour,  is  made  dimensionless  using  the  axial  dis¬ 
tance  between  the  nozzle  throat  section  and  the  exit  section. 
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Three  of  the  dimensionless  quantities  in  Eqs.  (43)  and  (44)  are 
defined  as  follows: 


(45) 

(46) 
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The  local  radius  of  the  Isentroplc  contour  has  been  made  dimen¬ 
sionless  using  the  radius  of  the  Isentroplc  contour  at  the 
nozzle  throat  section.  The  density  and  velocity  at  a  given 
section  of  the  isentroplc  contour  have  been  made  dimensionless 
based  on  the  exit  density  and  velocity,  respectively.  The 
length  Reynolds  number  is  defined  as: 


P  u  L 

Re  -  ~  C  -  (48) 

Lt  Uo 

Once  more  it  will  be  noticed  that  the  exit  density  and  velocity 
have  been  used  along  with  the  viscosity  evaluated  at  stagnation 
conditions.  The  dimensionless  axial  position  parameter,  x+,  is 
defined  as: 


AX 
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The  remaining  three 
as  follows: 
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quantities  in  Eqs.  (43)  and  (44)  are  defined 
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A  sine-arch  velocity  profile  in  the  boundary  layer  has  been 
employed  in  solving  the  integral  boundary  layer  equations. 

In  order  to  start  the  boundary  layer  calculation  at  the 
throat  section  of  the  nozzle,  it  is  assumed  that  the  boundary 
layer  growth  at  this  section  is  zero,  or: 


d6 

dx 


0 


(53) 


Using  this  condition  and  Eq.  (42),  it  is  possible  to  show  that 
the  boundary  layer  thickness  at  the  first  two  sections  along 
the  isentropic  contour  is  given  by: 


Once  the  boundary  layer  thickness  has  been  calculated  at  a 
given  section,  the  boundary  layer  displacement  thickness  can 
be  calculated  from  the  following  equation: 


n 


6+  L  (1  -  B  ) 
n  t  '  n7 


(55) 


The  actual  calculation  procedure  is  carried  out  as  follows. 
In  general  the  following  quantities  will  have  been  selected 
prior  to  the  design  of  the  detailed  nozzle  contour;  namely,  the 
exit  Mach  number,  an  isentropic  exponent  for  relating  the  exit 
conditions  to  the  stagnation  conditions,  the  stagnation  tem¬ 
perature,  and  either  a  stagnation  pressure  or  a  nozzle  exit 
pressure.  The  values  of  pressure  will  have  been  estimated 
through  the  use  of  charts  similar  to  Figs.  4  or  5  and  13  or  14. 
The  use  of  these  charts  will  have  given  preliminary  estimates 
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for  the  mean- free -path  at  the  nozzle  exit  section  and  will  have 
determined  vAiether  or  not  an  lsentroplc  core  will  exist  at  the 
nozzle  exit  section.  The  above  information  can  then  be  used  to 
calculate  the  following  quantities  at  the  nozzle  exit  plane; 
namely,  static  temperature,  density,  sonic  velocity,  actual  ve¬ 
locity  viscosity,  Reynolds  number  per  unit  length,  and  mean-free- 
path. 

Since  the  coordinates  for  most  lsentroplc  contours  are 
given  in  a  dimensionless  form,  it  is  necessary  to  pick  at 
least  one  important  dimension  for  the  lsentroplc  contour.  This 
can  be  done  through  the  use  of  the  approximate  boundary  layer 
analysis  and  usually  the  exit  diameter  for  the  lsentroplc  con¬ 
tour,  D  ,  will  be  chosen.  Using  this  dimension,  plus  the  area 
ratio  for  the  lsentroplc  contour,  it  is  impossible  to  determine 
the  throat  radius  and  the  quantity  Lfc.  At  this  point  all  of  the 
information  is  available  for  the  dimensionless  quantities  which 
appear  in  Eqs.  (43)  and  (44).  In  addition  to  the  actual  coor¬ 
dinates  for  the  lsentroplc  contour,  usually  values  of  Mach 
number  are  known  along  the  wall  of  the  lsentroplc  contour.  Since 
the  flow  field  for  the  lsentroplc  contour  corresponds  to  an 
lsentroplc  flow,  these  values  of  Mach  number  along  with  the 
stagnation  pressure  and  temperature  permit  the  calculation  of 
the  remaining  quantities  in  Eqs.  (43)  and  (44).  A  slightly 
different  approach  has  been  taken  for  the  nozzle  calculations 
\rtiich  will  be  presented.  Rather  than  employ  a  value  of  Mach 
number  given  along  the  wall  of  the  lsentroplc  contour,  the 
radial  coordinates  at  a  given  section  for  the  lsentroplc  contour 
have  been  used  to  calculate  the  area  ratio,  A/A*.  The  flow  has 
been  assumed  to  be  one-dimensional  and  these  area  ratios  have 
been  used  in  conjunction  with  one -dimensional  isentropic  flow 
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tables  to  calculate  the  required  density  and  velocities  In 
Eqs.  (43)  and  (44).  Reference  (52)  has  shown  that  this 
approximate  method  for  calculating  the  lsentroplc  flow  param¬ 
eters  does  not  Introduce  serious  errors  Into  the  calculations. 

After  the  boundary  layer  growth  along  the  nozzle  had  been 
calculated,  the  resulting  contour  was  plotted  on  large-size 
graph  paper  and  a  smoothing  process  was  employed  to  remove  any 
bumps  or  wavlness  which  had  been  Introduced  by  the  calculation 
procedure.  Only  the  smooth  nozzle  contour  Information  has  been 
sunmarlzed  In  this  report.  A  somewhat  more  detailed  procedure 
for  designing  nozzles  for  low-density  wind  tunnels  will  be  dis¬ 
cussed  later  In  this  section  after  the  calculated  results  for 
the  detailed  nozzle  designs  have  been  presented. 

Nozzle  Design  I 

The  first  detailed  nozzle  calculation  was  made  based  on 
the  assumption  that  nitrogen  would  be  used  as  the  working 
fluid  for  a  low-density  wind  tunnel.  The  stagnation  temperature 
was  also  maintained  at  a  fairly  low  value,  namely,  5500  R.  The 
use  of  nitrogen  at  this  temperature  led  to  a  set  of  conditions 
such  that  dissociation  of  the  nitrogen  was  not  an  important 
factor.  Therefore,  It  would  be  possible  to  avoid  non-equilibrium 

flow  effects  and  to  obtain  data  on  rarefied  gas  flow  phenomena 

without  worrying  about  the  complications  of  real-gas  effects.  An 

exit  Mach  number  of  10  was  selected,  and  an  exit  pressure  for  the 

nozzle  of  40  microns  was  used.  The  pressure  was  selected  based 
on  a  figure  similar  to  Fig.  4.  The  approximate  boundary  layer 
analysis  showed  that  there  should  be  some  lsentroplc  flow  region 
at  the  exit  of  the  nozzle.  These  calculations  are  given  below. 

At  the  time  this  nozzle  design  was  made,  only  one  lsentroplc  con- 


tour  corresponding  to  a  Mach  number  of  10  was  available  (50) . 
The  following  Information  was  specified: 

Nitrogen 

M  -  10 
e 

Y  -  1.4 

T  -  5500  R 
o 

pg  -  40  microns 

For  adiabatic  flow  of  a  perfect  gas,  the  exit  static  tempera¬ 
ture,  T  ,  can  be  calculated  as  follows: 
e 


1  +  0.2  (10) 2 


21 


T  .1500 
e  21 


261.9  R 


Once  the  static  temperature  is  known,  the  viscosity  and  sonic 
velocity  can  be  calculated  as  follows: 

0.073225  x  I0'5{7e  0.073225  x  10"5  Y 261.9 ' 

ue  "  201.6  "  201.6 

1  +  T  1  +  261.9 

€ 

U  ■  0.6695  x  10  5  lbm/ft-sec 
e 

a  -  fm  T  -  /  (32.174)  (1.4)  (1545.32)  (261.9) 

e  y  8o  Wo  e  r  (28) 

a  -  806.9  fps 
e  r 
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From  a  definition  of  Mach  number  the  exit  velocity,  Vfi, 
becomes : 

V  -  M  a  -  8069  fps 
e  e  e 

The  exit  density  can  be  calculated  from  Eq.  (7): 


P 

e 


W  p 
o  _ 

ZRT 


The  compressibility  factor  will  be  taken  as  1,  since  at  the 
temperature  involved  the  amount  of  dissociation  is  negligible 
especially  since  nitrogen  has  been  used.  Therefore  the  den¬ 
sity  becomes: 


(760,000) (1545.32) (261.9) 


-  7.704  x  10 


^lbm/ft"^ 


The  Reynolds  number  per  unit  length  can  be  calculated  from 
the  above  quantities  by  the  following  equation: 


Re/L 

-  773  inch'1 

The  mean-free-path  at  exit  conditions  is  calculated  from  Eq. 

(8): 

\  8  HP  Me  8  /  2(1.4)  *  10 

Ae"  5  ~  Ri7L“  773 

-  0.0195  inch 
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The  re suit 8  from  the  approximate  boundary  layer  analysis 

are  as  follows.  An  exit  diameter  for  the  isentropic  contour, 

D  ,  is  selected  as  A  Inches.  This  value  was  selected  based 
se 

on  the  results  shown  in  Fig.  12  since  Repge  will  be  about  3090 
using  the  calculated  Reynolds  number  per  unit  length.  This 
should  produce  a  value  of  6/r^  somewhat  less  than  0.9.  A  check 
on  the  isentropic  contour  (50)  will  also  show  that  0  **6°  so 
that  Fig.  12  could  be  used.  From  the  isentropic  contour  (50) 
the  quantity  Lt  becomes  38.28  inches.  From  Eq.  (19)  the  exit 
boundary  layer  thickness  is: 


6 


e 


C 


Ip  V  L 
feet 

V 


The  constant  C  in  Eq.  (19)  was  taken  from  Fig.  10  correspond¬ 
ing  to  the  Van  Driest  solution  for  adiabatic-wall  conditions. 
The  value  of  the  constant  C  is  31.8.  Putting  this  information 
into  the  above  equation,  the  boundary  layer  thickness  becomes: 


Also, 


6 

e 


the 


.  (31.8H38.28)  .  (31.8H  38.28)  ,  ?  Q6  lnche8 

J 773  x  38.28  J 29.600 

displacement  thickness  at  the  nozzle  exit  section  be¬ 


comes: 


6*  «  (1  -  B  )6 
e  e  e 

6*  -  (1  -  0.192)7.06  -  5.70  inches 
e 
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TABLE  III  -  CALCULATED  RESULTS  FOR  NOZZLE  DESIGN  I 
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The  actual  exit  radius  for  the  corrected  nozzle  contour  can  be 
calculated  from  Eq.  (34): 

r  -  r  +  6*  -  2  +  5.70  ■  7.70  Inches 
e  se 

The  radius  of  the  Isen tropic  contour  Is  determined  by  Eq.  (35) : 

r*  -  r  -  6  »  7.70  -  7.06  -  0.64  Inch 
se  e 

This  means  that  an  Isentroplc  flow  region  with  a  diameter  of 
approximately  1.28  Inches  will  be  available  for  test  purposes. 
It  will  be  seen  from  the  actual  nozzle  calculation  that  this 
estimate  was  over-optimistic. 

Table  III  summarizes  the  Important  information  for  the 
final  nozzle  calculation.  The  radius  of  the  isentroplc  con¬ 
tour,  the  displacement  thickness,  and  the  radius  of  the  actual 
contour  are  given.  Values  for  several  flow  parameters  are 
also  presented  including  the  Mach  number  distribution  for  the 
nozzle.  The  actual  boundary  layer  displacement  thickness  at 
the  exit  plane  of  the  nozzle  is  about  8.74  Inches.  Since  the 
radius  at  the  exit  section  of  the  isentroplc  contour  was 
taken  as  2  inches,  the  actual  exit  radius  for  the  final  nozzle 
is  10.74  inches.  Using  Eq.  (23)  to  convert  the  boundary  layer 
displacement  thickness  to  a  boundary  thickness,  yields  a  value 
of  the  boundary  layer  thickness  of  10.82  inches.  The  detailed 
nozzle  calculation  indicates  that  there  is  no  isentroplc  flow 
region  at  the  exit  plane  of  the  nozzle.  Therefore,  the  approx- 
iaate  boundary  layer  analysis  predicted  an  exit  diameter  for 
the  isentroplc  contour  tfilch  was  too  small.  It  is  interesting 
to  observe  that  the  boundary  layer  growth  at  Mach  numbers  close 
to  the  exit  Mach  number  is  extremely  rapid.  At  a  position  in 
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the  nozzle  located  halfway  between  the  nozzle  throat  and  the 
nozzle  exit  plane,  the  Mach  number  is  already  9.85,  and  the 
boundary  layer  displacement  thickness  is  5.3  inches.  This 
value  of  the  displacement  thickness  should  be  compared  with 
that  given  by  the  approximate  boundary  layer  analysis  above. 

In  the  downstream  half  of  the  divergent  section  of  the  nozzle, 
the  Mach  number  changes  from  9.85  to  10.068,  a  2  percent  change, 
while  the  boundary  layer  displacement  thickness  changes  from 
5.3  to  8.74  inches,  a  40  percent  change.  It  is  worthwhile  and 
important  to  determine  a  value  of  C  from  the  detailed  nozzle 
calculation.  This  can  be  done  using  the  equation  which  defines 
C: 

i  .  C 

Lt  "f.vA' 

The  value  of  6  from  the  detailed  boundary  layer  analysis  is 
used.  For  the  first  nozzle  design  C  is  48.8.  A  later  section 
is  devoted  to  further  comments  on  this  value  of  C,  which  is 
considerably  higher  than  that  given  by  the  Van  Driest  results 
for  adiabatic  flow  shown  in  Fig.  10. 

If  this  nozzle  contour  were  used  in  an  actual  low-density 
wind  tunnel,  under  the  conditions  assumed,  there  would  be  no 
region  of  uniform  isentropic  flow  in  which  to  place  a  model. 

Nozzle  Design  II 

The  second  detailed  nozzle  calculation  was  performed  using 
air  as  the  working  fluid.  The  stagnation  temperature  and 
pressure  were  6730  R  and  17.7  atm,  respectively.  The  exit  Mach 
number  for  the  nozzle  was  assumed  to  be  10.  Additional  isen¬ 
tropic  nozzle  contours  were  obtained  from  the  Sandberg-Serrell 
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Corporation  (51) .  A  nozzle  contour  corresponding  to  a  value 
of  1.2  was  used  for  the  second  nozzle  calculation.  Calculations 
pertaining  to  the  second  nozzle  contour  were  made  in  a  manner 
identical  to  those  for  the  first  nozzle.  The  results  are 
summarized  in  Tables  VI  and  VII. 

Table  IV  gives  the  nozzle  coordinates  for  the  second  con¬ 
tour.  The  boundary  layer  displacement  thickness  at  the  nozzle 
exit  is  14.54  inches.  The  corresponding  boundary  layer  thick¬ 
ness  is  19.50  inches.  Since  the  exit  radius  for  the  isentropic 
contour  is  only  2  inches,  there  is  no  isentropic  core  region 
for  this  second  contour.  This  situation  was  anticipated  by  the 
approximate  analysis  as  can  be  seen  in  Table  VII.  The  values 
for  the  displacement  thickness  from  the  detailed  and  approximate 
analyses  are  14.54  and  17.42  inches,  respectively.  For  the 
approximate  analysis,  a  value  of  C  of  25.9  was  obtained  from 
the  Howarth  result,  Eq.  (36),  using  a  value  of  v  of  1.2.  This 
value  of  C  predicted  values  for  the  boundary  layer  thicknesses 
which  were  too  large  in  comparison  to  the  values  from  the  de¬ 
tailed  analyses.  The  value  of  C  for  the  detailed  analyses  was 
only  21.6.  However,  the  results  were  sufficiently  close  enough 
to  each  other  so  that  the  use  of  the  approximate  analyses  can  be 
justified. 
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TABLE  IV  -  CALCULATED  RESULTS  FOR  NOZZLE  DESIGN  II 
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TABLE  IV  (continued) 
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Nozzle  Design  III 

The  third  detailed  nozzle  calculation  was  made  using  air  as 
a  working  fluid.  The  stagnation  temperature  and  pressure  were 
12500  R  and  89  atm,  respectively.  An  exit  Mach  number  of  10  was 
selected,  and  an  isentropic  contour  corresponding  to  a  value  of 
v  of  1.2  was  used  (51).  The  calculations  for  the  third  nozzle 
were  made  as  shown  for  the  first  nozzle  and  the  results  are 
given  in  Tables  VI  and  VII. 

Table  V  summarizes  the  nozzle  coordinates  for  this  third 
nozzle.  It  will  be  seen  that  the  exit  displacement  thickness 
and  the  boundary  layer  thickness  from  the  detailed  calculation 
are  9.47  and  12.72  inches,  respectively.  Since  the  exit  radius 
for  the  isentropic  contour  was  taken  as  2  inches,  once  more  the 
isentropic  flow  area  does  not  exist  at  the  exit  of  the  nozzle. 

For  this  third  contour  the  approximate  boundary  layer  analysis 
also  predicted  the  absence  of  the  isentropic  flow  region.  How¬ 
ever,  as  previously  encountered  in  the  second  nozzle  design, 
values  for  the  displacement  thickness  from  the  approximate  and 
the  detailed  calculations  are  very  close  to  each  other.  The 
value  of  the  displacement  thickness  from  the  approximate 
analysis  at  the  nozzle  exit  plane  was  11.27  inches  while  the 
detailed  calculation  gave  a  value  of  9.47  inches.  The  value  of 
the  constant  C  corresponding  to  the  detailed  calculation  was 
21.9  as  compared  to  the  Howarth  value  of  25.9.  Note  that  the 
values  of  C  from  the  detailed  analyses  of  the  second  and  third 
nozzles  are  approximately  equal. 

Summary 

Some  of  the  important  results  from  the  three  detailed  nozzle 
calculations  have  been  summarized  in  Tables  VI  and  VII.  The 
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general  operating  characteristics  for  the  low-density  wind 
tunnels  using  the  three  nozzles  are  given  in  Table  VI,  where 
the  important  test  section  conditions  have  been  presented 
along  with  some  of  the  tunnel  requirements  such  as  flow  rate, 
heating  requirements,  and  so  forth.  The  equations  and  dis¬ 
cussion  of  these  requirements  are  given  in  the  section  of 
this  report  entitled,  "Components  for  Low-Density  Wind  Tunnels". 
As  the  stagnation  temperature  was  increased,  at  a  constant 
value  of  Mach  number,  the  exit  velocity  increased  from  about 
8000  to  15000  ft/sec.  Since  a  different  value  of  Y  was  used 
for  two  of  three  nozzles,  the  static  pressure  requirements  do 
not  follow  a  regular  pattern.  In  general,  static  pressures 
between  7  and  40  microns  are  required.  The  mean  free  path  at 
test  section  conditions  varies  between  approximately  0.01  and 
0.313  inch. 

Table  VII  summarizes  the  boundary  layer  calculations  from 
the  approximate  and  from  the  detailed  calculations.  Consider 
the  first  detailed  nozzle  design  for  a  value  of  v  of  1.4.  It 
will  be  seen  that  the  use  of  the  Van  Driest  value  of  C  for 
adiabatic- flow  conditions  resulted  in  a  value  of  the  boundary 
layer  thickness  which  was  only  65  per  cent  of  that  obtained 
from  the  accurate  boundary  layer  analysis.  Instead  of  the  Van 
Driest  value  of  C  of  31.8,  a  value  of  C  of  48.8  would  be  re¬ 
quired  to  explain  the  value  of  boundary  layer  thickness  from 
the  detailed  analysis.  An  indication  that  such  a  large  differ¬ 
ence  might  be  expected  was  given  by  the  information  in  Table 
II,  where  values  of  C  from  existing  low-density  wind  tunnels 
were  compared  with  tne  values  from  the  Van  Driest  adiabatic- 
wall  solution.  Although  the  experimental  data  only  extend  up 
to  a  Mach  number  of  6,  at  this  Mach  number  the  experimental 
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value  was  about  1.5  times  larger.  The  values  for  C  from  Table 
II  and  from  the  first  nozzle  design  at  a  Mach  number  of  10  have 
been  plotted  in  Fig.  26,  together  with  values  from  the  Van 
Driest  adiabatic-wall  solutions  and  the  Howarth  solution  It 
will  be  seen  that  the  suggested  use  of  the  Van  Driest  solution 
is  fairly  accurate  at  Mach  numbers  of  about  2,  whereas  as  the  Mach 
number  is  gradually  increased  into  the  hypersonic  range,  the 
difference  between  the  Van  Driest  solution  and  the  average  curve 
through  the  existing  experimental  and  theoretical  data  for  a 
value  of  V  of  1.4  continues  to  increase  until  at  a  Mach  number 
of  10,  the  difference  is  about  60  percent.  It  will  be  seen 
that  the  predictions  from  the  Howarth  solution  for  a  value  of 
Y  of  1.4  are  in  much  better  agreement  with  the  mean  curve  for 
this  same  value  of  Y  than  those  from  the  Van  Driest  solution. 

Thus  the  mean  curve  is  available  for  values  of  C  when  y  is  1.4 
and:  the  Howarth  solution  may  be  used  as  a  close  approximation. 

The  values  of  C  from  the  second  and  third  detailed  nozzle 
designs  are  also  shown  in  Fig.  26.  Since  these  designs  corre¬ 
spond  to  a  value  of  Y  of  1.2  the  Howarth  solution  for  this 
value  has  also  been  shown.  The  closeness  of  the  results  at  a 
Mach  number  of  10  add  support  to  the  use  of  the  Howarth  solution 
for  obtaining  values  of  C  for  the  approximate  boundary  layer 
analysis  when  experimental  data  or  detailed  calculations  are 
absent. 
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TABLE  V  (continued) 


O^vO^MOOOvOincO  HO^OO 

«  oooooa>o>aNO>aN<^oooo 

COmrOCOrOMCNCNCMMCMCNN 


o 


cOOOsj’NrOino^vfi^  OOOVOCM 

n^vonHooNinstfONHo^ON 

oors^Nrs^)v£)voiovovo^^inin 

CsiCNlCSICNCSlCN4CSlCSlCNCNiCSlCN<NCMCNJ 


O 


OOOOOH^DHOOsfCMOSONHH 
iTlCOMOOOMTi^fOHOOaN 
<C  NfssSvOvOvOvOvOvO^OvO»T| 
CMCMCMCNCMCNCnJCNCsICNICSCNICsI 


o 


m^^HON^voo^NONcn 

HN^400in-sf 

mfnmNHvOMN^HOO^^ 

+  w 

a.  or^oorvr^vDvoiointnvtvt^ 


oNincncooooooooooooHoost 
in  csi  o>  in  ro  hoonoo  NNvomin 
lomvt^^^^^cncororoconco 
CNCSCSiCSICNCMCSCMCNCMCNJCMCNCNJCSJ 


O 


vomvo^ocsvoiosoP^asH^oovoro 

movCfOr-ioor^vomvJ'vJ’rocMcsicsi 


COCOCMCMCNiHHHHHHHHrlrl 


oovfoo^mN^NOinrnoo^ 
a\HN^invONOO(^OHH(M 
*fw  ^oNNNNrsrsNrsoooooooO 

2  C>  Os  Os  Os  Os  Os  Os  Os  Os  Os  Os  Os  ON 


o 


4  ro  o  in  rs  M  O^CSlOONf^^vON 
mi^CAOr-ico^<ttninv£)^or^r^r^ 
oooooocfNONONONO>ONON<?NONONcrNcr> 
On^OnOnON^OnO>CTn  Os  Os  Os  Os  Os  Os 


O 


ooNinsJ-HQOvtHNcooocMrN  oocnr^oNO^oinr^r^cnoocoinr^ 

vDNOO^OOHNNcnn^vj-  ^OOOOSONfO^^miTivOvONNN 


NNNNoOOOOOOOOOOOOOOOOO  OOCOCOOsOsOSOsOsOsOSOsOsOsOsOS 


c  O^OOoONr^vOvOvOvOvONNN 
•r4  CNfA^lAvONOO^OHCSm^ 


U  HHHHHHHHMWCSICMCM 


vOrsoOvOnON<AvOM^vO^OvOM 

QOCNvOO^NH^NO^vO^H^ 


csifOfO^^^jnminvDvovovDNrN 


* 

«o 


CNstvOOOO^vANO^inoOO 

r^cNcsiCNicooNcncn^r^o^o 

HNfoovOMoNvocnooomn 

vovor^aoao<^asOH<NCMro<t 


O  rH  iH  rH  rH  iH 


HHCSfOvO^NONCOOOOOMOON 

HCSCOOSH^ONOOvO^lOvOCO 

cncncococor^o^rHCMCMr^mvor^ 

so(Avoow«nNonvooooMin 


HCNcsMcnncoro^^^vj-ininin 


•  CMCOCOCM^OlA^QOlANfA^ 
C  OnOHN  HHONN^HNcnoO 
•H  CMvOa>CMinoOOfOvOONr-»<tvO 
>»/  \OvOvONNNOOOOQOOOO>ONO> 
<0 

o 


OOvOiAHOONvOCOflONGOONVOcno 

r0r-iCMr^inaNr^r-4O^orvmocMr-i 

vomcoor^cnoNino^oocMxooNCS 

OHCM(Am^<t»nv£)vO^NNSOO 


HrlHHHHrlrirlHHHHrlH 


•S 

K 


VOOOOCM^VOOOON'^'OOOO  oOvOsfCMOOOvO^NOflOO^NO 
mminminvONOsONOvON  r^ooasoHHCMco^-minvor^oooN 


74 


iHvOCMf^roosm^CSOONr^ 
flONN\c^inmininin^4 
inmtnmmmmmmmmm 
pa  CNCNCMCMCMCSiCNCNICMCMCsICSI 


o 


voHooincMcomcocNooNONoor^ 

<f^nrOcn(NCsi(MCNICSlHHHH 

cncocncncncncncncncncncncncnr  r  r 

<  CNCMCMCNCMCMCMCsICMCMCMCMCMCM 


O 


m  mo 

OOmnH^MTi^fnCJHHHOO 
+  CO  iHHfHHOOOOOOOOOOO 


CNirjvOOOOCSCO^vDNNNOOONO 

ooooooooonon^on^o^ono>o>o>o 
-f-  CO  onononononononononononononono 

G  On  ON  On  ON  On  ON  ON  ON  On  On  On  On  On  On  O 


S 


>3*  00  CM  VO 

cMinrs^HcnmvONooooooasON 

OOOOOOCJOONONOnONONONONONONONO 


ON  ON  ON  ON  ON  ON  On  On  ON  ON  On  ON  On  On  O 


G  ooomoNC\i<fmvO\o<tcMoomcNir^cnoNmomo  ^ 

1-1  V£)ONr-ICnvOOOOCN*>^NOONaNrHCO-4’V£)rvONHCNsd’  <t 


u 


r^r^OOOOOOC300NONaNONONONOOOOOOrHr-lr-<  1—1 


.  inHr^csioNi-<r^cnor-ior>*mmcMo<i*oiH<i*<M  cn 

C  r^HmmoNinaNrocommoNCN40Nvoojr^r-icocor-i  no 

•W  ONOmOO>frCMr^r^sO<fCNO<Mr-ir^or)OOvtON<J’ON  vo 

w  r^o<M>tNOONOCS<fvooNOiHro<tvor^ONOCNm 


o  invDvoovovor^1'NrsNNoooooocooocooo(j\ONON  on 


T) 

0) 

'd 

G 


o 

G 


o 

o 


.  ONCOCsiOOCOoOCO.HCOONCMCMrHOOmHmr^OOCO  oo 

G  vOOl— lONNOOm^COO^CMvOONrlCMfOCOCOCOCO  cn 

•h  <tr>ONOCN<ttnor>.oooooNONO>ooooooo  o 

w  OOOOOOONONONONONONONONONONONOOOOOOO  O 

CO . 

^  r— H  rH  i— I  r-H  nH  r*H  i—H  rH  r~H  p-H  fH  r-H  pH  (Nj  CN  (N  ON  CM  CM  CM  CN 


W 

* 

f5 


C  OOvO^CNOOOs£)4NOOO'0<t<SIOOO'i)^CNOOO  CM 


ONOH<Mnco<f»nvorsNoooNOHHCNjr)^inin  vo 
X  r- ICMCMCMCMCMCMCMOsICNCMCMCMcncncncOCOcncncn  cn 


75 


vO 


.0) 


oc5 


co  *3*  vo 

r-  <f  o 

rH 


25 

o 

H 

H 

O 

W  W 
c/5  Z 

o 

H  H 

CO  CO 

w  w 

H  O 

Q  W 

i  ^ 

Csi 

B  § 

1  § 
a  m 
H  < 
P  H 
&  W 

a : 

fa  S 

5  co 
P  5s 

H  O 


>  I 

w  u 


m 

H 


uo 

X  |  o> 

CO 

rH 

0  H 

H 

CO 

c  o 

CO 

rH 

•H 

• 

• 

O 

O 

O 

o 

o 

C0  1  O 

o 

o 

a  I  Is* 

CM 

CO 

4-1  O 

uo 

1  00 

rH 

rH 

rH 

1  CM 

CM 

vO 

pet  vo 

rH 

co 

CM 

vO 

r-i 

(0 

C 

o 

0)  Vi 

0*  O 

iH 

a 


s  i 


H°  Pei 


X°  fl 


o 


vO 

r*-  oo 
co 


o  o  o 


o 

o 

o 

o 

o 

CO 

m 

m 

in 

vo 

CM 

H 

CM 

O 

• 

• 

« 

CM 

a> 

00 

CM 

CM 

• 

• 

• 

H 

|H 

rH 

co 


o 

o> 

co 

CO 

M 

Q> 

4J 


I 


o 

Q> 

CO 

J 


J3 
<u  o 

°  5 


^3 

a 

5 


o  o  o 

o  o  o 

o  o  o 

©  r**  cm 

CM  CM  CO 


CM  ^  o 

•  •  t 

ON  H  O 


o 

rH 

o 

<f 

vO 

o 

m 

O 

CM 

o 

o 

O 

o 

o 

O 

• 

• 

o 

o 

o 

CO 

CO 

vf 

in 

• 

• 

• 

o 

vO 

H 

rH 

rH 

H 

>4 

CM 

CM 

CM 

CM 

rH 

O 

O 

• 

• 

• 

O 

O 

O 

CO 

CO 

a 


c 

01 

00  V  V 

O  *rl  *iM 

u  <  < 

4J 


<u  pd 

co  a 

Q  C 


<4-  ^  ^ 


N 

N 

o 

25 


H  M 


V 

i— I 

N 

N 

O 

2 


H  M 


M  H 


76 


o 

cn 

CO 

rH 

CM 

CM 

00 

vO 

ON 

00 

C\ 

on 

• 

• 

• 

• 

• 

• 

00 

H 

1— 1 

H 

uo 

uo 

Mt 

CM 

CM 

CO 

CM 

CM 

CD 
•  CO 

Q 


43 

0 

c 


VO  CM  O 

i-i  as  m 


o 

i 


m  cm 

i  i 


00  VO  O 

CM  00  P** 

•  •  • 
H  P^  CO 
+  •  » 


5 

CO 

CO 

z 

a 

pi 

< 

£ 

2 

& 

< 

Q 


a 

2 

Q 

W 


•g 

c 


*£3 
a)  u 
a  c 


CM 

o 

CM 

vO 

in 

CM 

00 

m 

P*- 

O 

CO 

rH 

• 

• 

• 

0 

• 

o 

ON 

CM 

CO 

m 

*rH 

rH 

rH 

CM 

rH 

00 

00 

Mt 

O 

o 

ON 

vfr 

00 

m 

# 

• 

• 

• 

• 

rH 

CO 

CM 

m 

00 

NO 

CM 

CO 

CM 

rH 

CO 

CM 

«  e 

z  o 
o 

M  Q 
H  Z 
a  < 
w 

w  g 

H  < 


Q  g 

.  § 
a, 

H  Pm 
H  < 
> 


w 

a 

?5 


42 

a 

c 


CO 

s 


< 

3 

o 

§ 

w 

ri 

tS 

W 

P 


<D 

(0 


43 

O 

<3 


0) 

r-l 

N 

N 

O 

Z 


o 

CM 

m 

p^ 

« 

♦ 

• 

• 

• 

00 

rH 

ON 

m 

rH 

<t  ^ 


CO 

§ 


i 

B 


H  >< 
H  9 

H  Pi 
Pm 
Pm 
< 


CM 


<fr  Mt  <t 


h  cm  co 


77 


From  Fig.  10,  Van  Driest  Adiabatic-Wall  Values 
From  Howarth,  Eq.  (36) 

Values  from  (6/l.)  -  C/  ip  VI./u  '  using  values  of  6  from  detailed  calculation 


SUGGESTED  DESIGN  PROCEDURE  FOR  LOW-DENSITY  WIND  TUNNELS 


In  view  of  the  results  of  the  approximate  and  detailed  ana¬ 
lysis,  it  is  possible  to  identify  a  design  procedure  for  low- 
density  wind  tunnels.  It  is  assumed  that  at  the  start  of  the 
design  two  quantities  are  specified: 

Tq  -  stagnation  temperature 

M  -  nozzle  exit  Mach  number 
e 

Equation  (27)  should  be  used  to  determine  values  of  Repse  as  a 

function  of  */ r  . 

e 


Starting  with  a  value  of  6/re  of  unity  and  selecting  subsequent 

values  based  on  a  fine  increment,  i.e.,  0.98,  0.96,  etc.,  a 

series  of  values  of  Ren  can  be  obtained.  There  are  two  reasons 

Dse 

for  using  a  fine  increment.  First,  the  size  of  the  nozzle  grows 
very  rapidly  as  the  increment  is  increased  and  the  sizes,  there¬ 
fore,  may  become  ridiculous.  Secondly,  the  use  of  a  fine  incre¬ 
ment  will  yield  a  sufficient  amount  of  data  on  the  expected  is- 

entropic  flow  diameter,  D'  ,  so  that  a  decision  can  be  made  *rtien 

se 

this  region  has  achieved  a  size  so  that  uncertainties  in  the 
approximate  analysis  will  not  cause  it  to  vanish  when  the  detail¬ 
ed  calculation  is  made.  The  remaining  quantities  in  Eq.  (27) 
are  to  be  selected  as  follows: 

1.  I sen tropic  Exponent,  Y  -  The  value  of  Y  in  Eq. 

(27)  is  used  to  relate  thermodynamic  properties  at 
the  Mach  number  Mg  to  stagnation  conditions.  Thus 
from  the  various  plots  given  in  this  report  on  nozzle 
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expansion  processes  it  is  obvious  that  the  values 
of  v  will  vary  between  1.1  and  1.67.  Limiting  the 
discussion  to  expansions  to  high  Mach  numbers,  i.e. 
10  or  greater,  it  appears  that  for  tunnels  using 
air  at  stagnation  pressures  in  the  range  from  1  to 
100  atmospheres  and  stagnation  temperatures  around 
5000K,  a  value  of  Y  of  about  1.4  should  be  used  in 
Eq.  (27).  For  a  diatomic  gas,  the  ease  with  *rtiich 
the  gas  can  be  dissociated  will  exert  a  strong  in¬ 
fluence  on  the  value  of  Y.  At  a  given  pressure 
and  temperature,  oxygen  may  be  in  a  monatomic  form 
while  nitrogen  may  be  mostly  in  a  diatomic  form. 

In  the  general  case  it  will  be  necessary  to  make 
calculations  on  the  real  expansion  processes  and 
from  these  deduce  a  value  of  Y. 

2.  B  -  Equation  (24)  should  be  used  with  the 

known  values  of  v  and  M  . 

e 

3.  6  -  A  value  of  6°  seems  representative  for 
most  isentropic  contours. 

4.  A/A*  -  Values  of  the  area  ratio  follow  im¬ 
mediately  from  the  known  values  of  Y  and  Mfi. 

5.  C  -  The  average  curve  given  in  Fig.  26  should 
be  used  for  low-density  wind  tunnels  using  air  in 
the  range  of  stagnation  pressures  and  temperatures 
mentioned.  For  other  gases,  detailed  nozzle  de¬ 
signs  must  be  done  before  the  values  of  C  are 
available.  Based  on  the  results  shown  in  Fig.  26, 
it  appears  that  the  Howarth  solution,  Eq.  (36), 
will  predict  reasonable  values  for  C. 
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Values  of  the  Reynolds  number  per  unit  length  should  now 
be  selected  and  the  actual  diameters  and  boundary  layer  thick¬ 
nesses  should  be  calculated  from  the  following  results: 


se 


ReDSe/(Ee/L) 


(Re/L) 


(* 


'6/r 


-  1)  +  1 


1 


(33) 


(30) 


6  -  6  (1  -  B  ) 

e  e  e 


D  -  D  +  26 
e  se  e 


(23) 

(34) 


D'  -  D  -  26 
e  e 
se 


(35) 


The  results  should  then  be  Inspected  in  order  to  obtain  the  low¬ 
est  possible  value  of  Re/L.  This  minimum  value  will  yield  the 
largest  mean-free-path  and  will  allow  the  use  of  the  largest 
model  at  a  given  Knudsen  number.  In  general,  the  largest  nozzle 
size  will  be  selected.  The  isentropic  flow  region  must  be  suf¬ 
ficiently  large  and  the  nozzle  must  be  compatible  with  facility 
supply  systems  and  budgets. 

The  test  section  pressure  and  mean-free-path  and  the  stag¬ 
nation  pressure  can  now  be  calculated  from  the  following  results 


s..  s  p  /i ITT 

L  7*1. 


(ii) 
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The  facility  requirements  can  be  calculated  from  Eqs.  (57),  (58), 
and  (61). 

Detailed  calculations  for  the  nozzle  can  now  be  made.  Due 
to  the  variation  of  the  isentropic  exponent  in  the  frozen-flow 
expansion  from  the  nozzle  throat,  an  isentropic  contour  for  the 
nozzle  should  be  calculated  taking  into  account  the  variation  of 
thermodynamic  properties.  Once  this  contour  has  been  determined, 
it  is  possible  to  calculate  the  actual  nozzle  expansion  process 
for  this  particular  contour  using  present  data  on  the  dissociation 
and  recombination  rates.  This  check  on  the  isentropic  contour 
design  will  reveal  the  point  in  the  expansion  process  at  which 
freezing  actually  occurs  and  will  indicate  whether  a  new  isen¬ 
tropic  contour  should  be  designed.  The  boundary  layer  corrections 
can  then  be  calculated  using  either  the  method  of  this  report  or 
based  on  the  results  of  several  papers  which  appeared  in  the 
literature  after  this  study  had  been  completed;  namely.  Refs. 

(53),  (54),  and  (55).  In  the  event  that  the  approximate  analysis 
had  been  overoptimistic,  in  the  prediction  of  the  boundary  layer 
thicknesses,  it  may  be  necessary  to  repeat  the  isentropic  contour 
and  boundary  layer  calculations  until  a  sufficiently  large  core 
of  isentropic  uniform  flow  is  produced  at  the  nozzle  exit  section. 
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COMPONENTS  FOR  LOW-DENSITY  WIND  TUNNELS 


There  are  three  important  components  for  low-density  wind 
tunnels  which  have  not  been  considered  thus  far  in  this  report. 
These  components  are  the  heater,  the  precooler,  and  the  low- 
pressure  vacuum  source  for  the  wind  tunnel.  Although  time  was 
not  available  for  a  complete  study  of  these  components,  some 
comments  can  be  made  based  on  some  rather  simple  and  accurate 
calculations. 

Mass  Flow  Rate  Requirements 

The  calculations  of  nozzle  expansion  processes  for  air 
over  the  stagnation  temperature  and  pressure  ranges  of  interest 
indicate  that  the  expansion  process  from  the  stagnation  state 
to  the  nozzle  throat  will  be  an  equilibrium- flow  process.  The 
isentropic  exponent  for  these  processes  will  have  values  of 
about  1.2  or  less.  The  compressibility  factor  will  be  about 
1.25  for  stagnation  pressures  between  1  and  100  atmospheres  if 
the  stagnation  temperatures  do  not  exceed  about  6500  K.  (See 
Fig.  1,  Ref.  5)  The  choked  flow  rate  can  be  calculated  using 
the  usual  expression  for  one -dimensional  flow  of  a  perfect  gas: 


w 

A* 


14.7  p 
_ |_o 

y*r 


where 


w  -  lbm/ sec 
A*  -  in2 

gQ  -  32.174  lbm  ft/lbf  sec2 

R  -  1545.32  ft  lbf/lbm-mole  R 

pQ  -  atmospheres 

T  -  R 
o 


(57) 
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The  differences  between  the  choked  flow  rates  from  the  equiiib- 
r I um-air  calculations  and  those  from  Eq.  (57)  using  "appro¬ 
priate"  values  for  y  and  ZQ  will  not  exceed  a  few  percent. 

Choked  flow  rates  were  calculated  for  values  of  Y  of  1.1,  1.2, 

1.3  and  1.4.  Values  of  T  and  p  were  selected  and  then  the 

o  o 

corresponding  value  of  Zq  was  obtained  from  Reference  (56)  . 

The  results  are  shown  in  Fig.  27  for  Y  -  1.2.  It  was 
found  that  for  other  values  of  y,  the  choked  flow  rates  were 
about  6%  higher  for  Y  -  1.4  and  3%  lower  for  Y  -  1.1.  The 
mass  flow  rate  per  unit  throat  area  has  been  given  for  an  ob¬ 
vious  reason.  At  a  stagnation  temperature  of  3000R  the  mass 

2 

flow  rate  varies  between  0.135  and  135  lbm/ sec-in  over  the 
stagnation  pressure  range  from  1  to  1000  atmospheres.  At  a 
higher  stagnation  temperature  of  14000R  the  flow  rates  are 
lower  by  a  factor  of  about  two.  At  a  constant  value  of  stag¬ 
nation  pressure,  the  mass  flow  rate  should  not  depend  on 
1/2 

Tq  since  Zq  is  also  a  function  of  stagnation  temperature. 

For  a  similar  reason  the  mass  flow  rate  should  not  necessarily 

be  a  linear  function  of  pQ  although  for  the  range  of  variables 

shown  in  Fig.  27,  the  predicted  linear  behavior  is  observed  at 

all  but  highest  values  of  stagnation  temperature. 

In  order  to  use  the  results  shown  in  Fig.  27  specific 

operating  conditions  for  the  wind  tunnel  must  be  selected. 

From  Fig.  13  the  exit  diameter  for  the  isentropic  contour  is 

1.36  inches  if  an  exit  Mach  number  of  10  is  required.  Since 

a  value  of  y  of  1.4  was  used  to  relate  the  nozzle  exit  state  to 

the  stagnation  conditions,  the  area  ratio  of  the  nozzle,  A /A*,  is 

2 

536,  D*  is  0.0587  inch  and  the  throat  area  is  0.00271  inch  .  This 
diameter  corresponds  to  the  minimum  exit  diameter  of  the  isen¬ 
tropic  contour  or  the  closing  of  the  boundary  layers  at  the 
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exit  of  the  nozzle,  l.e.,  6/r  ■  1.  Since  It  was  based  on  the 

Van  Driest  values  for  the  boundary  layer  growth  constant,  C, 

some  adjustment  should  be  made  for  the  differences  in  C  shown 

In  Fig.  26.  From  Eqs.  (28)  and  (33)  It  will  be  seen  that  the 

2 

diameter  depends  on  C  .  At  a  Mach  number  of  10  the  Van  Driest 

value  for  C  is  about  half  the  recommended  value  shown  in  Fig. 

26.  The  throat  diameter  should  be  increased  by  a  factor  of  4 

and  the  throat  area  should  be  increased  by  a  factor  of  16. 

The  new  values  of  these  quantities  then  become  0.235  inch  and 
2 

0.0434  inch  .  These  dimensions  must  now  be  related  to  other 

tunnel  conditions.  A  value  of  the  Reynolds  number  per  unit 
3-1 

length  of  10  inch  from  Fig.  13  must  be  retained.  For  a 
stagnation  temperature  of  12,000R  (Fig.  5)  the  static  pressure 
at  the  nozzle  exit  is  about  170  microns.  This  value  was  fixed 
by  the  specified  values  of  M,  Re/L  and  Y.  From  Fig.  6  the 
stagnation  pressure  must  be  9.35  atmospheres.  Finally,  the 
values  of  p  and  T  fix  the  mass  flow  rate  in  Fig.  27  at  about 

o  A  O 

0.55  lbm/sec-in  .  The  required  mass  flow  rate  is  0.0238  lbm/sec. 
Requirements  for  Heater 

The  minimum  power  requirements  for  the  heater  can  be  cal¬ 
culated  from  the  steady-flow  energy  equation: 

A*  -  <?*>  <ho  *  V  <58) 

The  mass  flow  rates  per  unit  area  were  taken  from  Fig.  27.  The 
values  of  h  were  taken  from  Reference  (56) .  The  term  h  was 

O  O 

evaluated  at  70F.  The  results  are  shown  in  Fig.  28  where  the 

2 

heater  power  is  given  as  kw/in  . 
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At  the  highest  stagnation  temperatures  and 

2 

pressures,  the  heater  power  is  about  100,000  kw/in  .  Since 

2 

throat  areas  between  0.01  and  0.1  in  are  to  be  expected,  the 
heater  power  requirements  are  probably  in  the  range  from  1000 
to  10,000  kilowatts.  For  the  example  given  in  the  first 
section,  the  stagnation  conditions  of  9.35  atmospheres  and 

12,000  R  require  a  heater  power  per  square  inch  of  about  3000 

2  2 
kw/in  and  the  throat  area  of  0.0434  in  results  in  a  heater 

power  of  about  130  kilowatts. 

Requirements  for  Precooler 

The  exact  requirements  for  the  precooler  will  be  deter¬ 
mined  by  the  specific  type  of  exhaust  system  which  is  employed. 
The  use  of  a  cryogenic  pumping  system  does  impose  some  obvious 
restrictions.  Presently,  the  largest  cryogenic  unit  being 

built  has  a  capacity  of  1300  watts.  The  lowest  heater  require- 

2 

ments  (for  the  use  of  air)  are  about  100  kw/in  (Fig.  29)  and 

2 

for  throat  areas  between  0.01  and  0.10  in  ,  the  heater  power 
inputs  become  1000  to  10,000  watts.  These  results  indicate  that 
the  precooler  for  a  low-density  wind  tunnel  using  a  single 
cryogenic  pumping  unit  must  have  a  capacity  equal  to  the  heater 
power  input  and  any  compressor  power  which  is  involved.  Cal¬ 
culations  were  made  for  air  which  showed  that  the  compressor 
power  was  several  orders  of  magnitude  smaller  than  the  heater 
power.  In  view  of  the  heater  requirements  and  the  capacity  of 
future  cryogenic  pumping  units,  it  appears  that  the  precooler 
must  always  have  the  capacity  just  mentioned.  The  design  of 
the  precooler  must  also  be  such  that  the  pressure  drop  is  made 
compatible  with  the  cryogenic  unit. 
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If  a  steam-ejector  system  is  used,  then  the  specifications 
for  the  precooler  are  indeterminate  until  the  characteristics  of 
the  ejector  system  are  known.  The  precooler  could  be  a  surface- 
to-gas  heat  exchanger  whose  capacity  was  adjusted  to  yield  sat¬ 
isfactory  operation  of  the  ejector.  The  operating  point  would 
be  determined  when  enough  heat  was  removed  from  the  gas  stream 
so  that  the  ejector  performance  was  still  acceptable  and  when 
the  precooler  pressure  drop  was  matched  with  the  wind  tunnel  and 
ejector  characteristics.  The  design  of  a  precooler  with  a 
sufficiently  low  pressure  drop  may  be  difficult  or  impossible  to 
achieve.  The  precooler  might  also  be  a  spray- tube  cooler  in 
which  cooling  water  is  introduced  into  the  gas  scream.  In  this 
case  a  suitable  operating  point  is  achieved  when  the  total  flow 
rate  of  gas  and  coolant  is  matched  to  the  ejector  characteristics 
and  any  pressure  drop  restrictions  have  been  met. 

Requirements  for  Exhaust  System 

Some  remarks  have  already  been  made  on  the  exhaust  systems. 
The  cryogenic  system  cannot  tolerate  any  heat  load  other  than 
that  required  to  condense  the  gas.  The  maximum  capacity  of  a 
single  unit  is  about  1300  watts.  Since  the  enthalpy  change  for 
air  between  the  saturated  vapor  and  saturated  solid  states  at 
low  temperatures  is  about  9  Btu/lbm,  the  mass  flow  rate  capacity 
of  a  single  unit,  with  no  additional  heat  load,  is  about  0.14 
lbm/sec.  This  number,  together  with  the  information  given  in 
Fig.  27,  then  define  the  allowable  throat  sizes  for  low-density 
wind  tunnels  using  this  single  cryogenic  unit. 

A  final  quantity  of  interest  is  the  volume  flow  rate  at 
nozzle  exit  density.  This  quantity  is: 
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(59) 


S 
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Using  Eq.  (57)  the  following  expression  is  obtained: 


(60) 


Z  .  ~  •  fT*  (61) 

e  °e  1  o 

For  the  calculations  the  following  values  were  used: 

v  -  1.2 

2  -  1.1 
o 

Z  -  1.0 

e 

The  volume  flow  rate  depends  on  the  stagnation  temperature  and 
the  exit  Mach  number  through  the  density  ratio.  A  value  of  v 
of  1.4  was  used  to  determine  the  density  ratio  since  this  value 
is  more  characteristic  of  the  overall  nozzle  expansion  process. 
The  value  of  V  of  1.2  was  used  for  the  process  between  the  stag¬ 
nation  state  and  the  nozzle  throat.  The  results  are  shown  in 
Fig.  29.  At  the  higher  Mach  numbers  the  volume  flow  rates 
approach  one  million  liters  per  second  per  square  inch.  Whether 
viewed  from  a  mass  flow  rate  or  a  volume  flow  rate  basis ,  the 
requirements  for  a  steam-ejector  system  are  enormous. 
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CONCLUSIONS 


The  most  critical  component  in  future  low-density  wind 
tunnels  which  are  designed  to  operate  at  high  stagnation  tem¬ 
peratures  is  the  nozzle.  These  nozzles  will  be  operated  with 
stagnation  temperatures  between  3000  and  10000  K  and  with  stag  - 
nation  pressures  between  0.1  and  100  atmospheres.  Due  to  this 
combination  of  pressures  and  temperatures,  the  actual  expansion 
process  in  the  nozzle  will  consist  of  three  regions:  the  ex¬ 
pansion  process  from  the  stagnation  state  to  some  section  down¬ 
stream  of  the  throat  will  be  an  equilibrium- flow  process  in 
which  the  chemical  reaction  rates  may  be  considered  as  infinite. 
For  some  distance  downstream  of  the  throat  a  second  region  will 
be  encountered  in  which  the  rates  of  dissociation  and  recom¬ 
bination  are  important.  This  will  be  a  nonequilibrium- flow 
region  in  which  the  chemical  kinetic  processes  must  be  consid¬ 
ered.  The  last  stage  of  the  nozzle  expansion  process  will  be  a 
frozen  flow  process  in  which  the  composition  of  the  gas  mixture 
will  be  constant  but  in  which  the  thermodynamic  properties  may 
be  either  variable  or  constant.  It  will  be  found  that  a  partic¬ 
ular  value  of  the  isentropic  exponent  will  be  required  to 
relate  the  nozzle  exit  section  conditions  to  the  nozzle  stag¬ 
nation  conditions,  while  an  average  or  varying  value  of  the  is¬ 
entropic  exponent  will  be  required  to  describe  the  actual  nozzle 
expansion  process. 

An  approximate  boundary  layer  analysis  has  been  developed 
tAiich  can  be  used  to  predict  feasible  operating  ranges  for  low- 
density  wind  tunnels.  The  results  of  this  analysis  lead  to  in¬ 
formation  on  the  nozzle  stagnation  pressure,  temperature,  and 
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nozzle  sizes  for  Isentroplc  contours  which  will  produce  various 
sizes  of  uniform  isentroplc  flow  at  the  nozzle  exit  section. 

In  order  to  utilize  the  approximate  boundary  layer  analysis,  it 
is  necessary  to  have  a  knowledge  of  the  boundary  layer  velocity 
profile,  a  value  of  ]f  to  describe  the  expansion  process  between 
the  nozzle  stagnation  state  and  the  nozzle  exit  section,  and  a 
value  of  boundary  layer  growth  constant,  C.  The  sine-arch 
velocity  profile  has  been  used  previously  in  the  design  of  low- 
density  wind  tunnel  nozzles  and  was  used  in  this  report.  The 
actual  nozzle  expansion  processes  which  will  be  encountered  in 
low-density  wind  tunnels  are  such  that  values  of  the  isentroplc 
exponent  between  1.1  and  1.67  are  to  be  expected.  Values  of 
the  boundary  layer  growth  constant,  C,  have  been  obtained  from 
existing  experimental  data  for  low-density  wind  tunnels  up  to 
a  Mach  number  of  about  6,  while  beyond  this  Mach  number  an 
extrapolated  curve  is  available  (Fig.  26) .  Values  of  C  at  a 
Mach  number  of  10  were  determined  by  detailed  calculations  dur¬ 
ing  this  study. 

Based  on  the  approximate  boundary  layer  analysis,  it  is 
possible  to  determine  the  required  static  pressure  at  the  exit 
of  the  nozzle  and  hence  the  stagnation  pressure  for  the  wind 
tunnel.  The  mean- free -path  at  the  exit  of  the  nozzle  can  also 
be  calculated,  and  this  information  can  be  used  to  determine 
the  allowable  ranges  of  model  size  and  Knudsen  number  which  can 
be  produced  by  the  tunnel.  It  is  also  possible  to  calculate 
the  power  requirements,  mass  flow  rates,  and  exhaust  system  re¬ 
quirements  for  the  tunnel. 

After  the  approximate  analyses  have  been  completed,  a  de¬ 
tailed  nozzle  design  must  be  made  in  which  an  isentroplc  contour 
is  generated  using  the  variable  thermodynamic  properties  which 
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characterize  the  actual  nozzle  expansion  process.  The  boundary 
layer  correction  calculations  can  then  be  made  based  on  the 
analysis  given  in  this  report  or  upon  publications  which  were 
made  after  the  completion  of  this  study. 

One  of  the  most  important  conclusions  of  this  study  is  re¬ 
lated  to  the  future  applications  of  low-density  wind  tunnels. 

It  does  not  appear  that  it  will  ever  be  possible  to  utilize 
"large"  models  in  low-density  wind  tunnels  under  conditions 
such  that  Knudsen  numbers  corresponding  to  free -molecular- flow 
conditions  can  be  generated.  Such  conditions  are  deemed  as 
being  impossible  due  to  the  huge  size  of  a  low-density  wind 
tunnel  which  would  be  required.  A  large  model  is  considered 
to  be  one  whose  characteristic  dimension  is  approximately  1 
inch.  The  future  of  low-density  wind  tunnels  should  be  devoted 
to  testing  the  largest  possible  models  which  can  be  used  under 
conditions  in  the  slip  and  transition  flow  regimes.  The  use  of 
large  models  will  permit  the  highest  possible  accuracy  for 
measurements  in  these  flow  regimes.  Such  accuracy  is  required 
in  the  development  of  the  theoretical  bases  for  these  flow 
regimes.  The  main  effort  on  research  in  the  free-molecular- 
flow  regime  should  be  devoted  to  basic  studies  of  the  surface- 
particle  interaction  phenomena.  Once  more  a  large  amount  of 
highly  accurate  data  is  required  for  the  surface-particle  in¬ 
teraction  phenomena  in  order  to  improve  the  theoretical  bases 
for  these  phenomena  which  are  presently  in  poor  condition. 

This  study  has  also  pointed  out  the  need  for  the  use  of  gases 
other  than  air  in  high- stagnation  hypervelocity  low-density 
wind  tunnels.  An  understanding  of  the  basic  processes  which 
are  Involved  between  a  rarefied  gas  stream  and  a  model  in  the 
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slip,  transition  and  free-molecular- flow  regions  is  difficult 
enough  to  obtain  without  the  complications  which  arise  from 
effects  such  as  a  nonequilibrium  flow  condition.  The  use  of 
gases  other  than  air  or  other  dissociating  gases  would  allow 
a  research  effort  to  separate  out  those  effects  which  are 
peculiar  to  rarefied  gas  dynamics  and  those  effects  which  are 
peculiar  to  high- temperature  physical  phenomena. 
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RECOMMENDATIONS 


Based  on  the  results  of  this  study,  it  is  possible  to  make 
the  following  recommendations  for  future  work  on  the  design  of 
low-density  wind  tunnels. 

Additional  calculations  should  be  made  on  the  actual  nozzle 
processes  which  will  be  encountered  in  low-density  wind  tunnels, 
since  present  calculations  of  this  type  exist  only  down  to  stag¬ 
nation  pressures  of  100  atmospheres.  The  calculations  should 
be  extended  to  stagnation  pressures  in  the  range  from  0.1  to 
100  atm.  The  stagnation  temperature  range  should  be  between 
3000  and  10000  K.  A  survey  of  available  data  on  the  dissociation 
and  recombination  rates  for  air  and  for  the  gases  of  which  air 
is  comprised  should  be  made.  Such  a  survey  is  required,  since 
considerable  effort  is  being  made  at  the  present  time  to  deter¬ 
mine  these  rates,  and  new  data  are  becoming  available  at  a  high 
rate.  In  addition  to  air,  other  gases  should  be  investigated; 
particularly  those  which  will  not  show  large  effects  of  non¬ 
equilibrium-flow  phenomena. 

Isentroplc  nozzle  contours  for  low-density  wind  tunnels 
should  be  designed  based  on  a  variable  isentropic  exponent. 

Some  procedures  already  exist  to  account  for  this  effect;  how¬ 
ever,  they  have  not  been  developed  to  the  point  where  they  are 
readily  usable  for  rapid  calculations.  Some  consideration 
should  also  be  given  to  the  design  of  nozzle  contours  when 
non-equilibrium-flow  conditions  are  encountered. 

A  review  should  be  made  of  accurate  boundary  layer  correc¬ 
tion  techniques  for  use  in  high- temperature,  low-density  nozzles. 
The  method  of  Reference  10  was  used  in  this  report;  however, 
additional  work  (53),  (54),  (55),  came  out  after  the  completion 
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of  this  study  and  these  new  methods  should  be  investigated.  Com¬ 
parisons  should  be  made  between  the  results  predicted  by  the  dif¬ 
ferent  methods  apd  some  conclusions  drawn  as  to  whether  the  state 
of  the  art  is  now  acceptable  or  \rtiether  further  work  is  required. 
In  relation  to  the  accurate  boundary  layer  correction  techniques, 
it  should  be  mentioned  that  a  study  should  also  be  made  of  the 
integral  boundary  layer  equations  when  dissociation,  ionization, 
and  so  forth,  are  present.  Although  a  solution  of  the  boundary 
layer  equations  using  the  integral  method  are  highly  developed 
for  incompressible  and  low-speed  flows,  additional  work  is  re¬ 
quired  when  dissociation  effects,  etc.,  are  present.  Existing 
data  and  theoretical  predictions  for  boundary  layer  separation 
in  low-density  wind  tunnel  nozzles  should  be  incorporated  into 
the  design  of  nozzles. 

The  approximate  boundary  layer  analysis  which  was  developed 
in  this  report  is  an  extremely  useful  tool  for  determining  fea¬ 
sible  operating  ranges  for  low-density  wind  tunnels.  It  is  also 
useful  in  preliminary  nozzle  design  estimates.  However,  it  is 
obviously  lacking  one  important  ingredient;  namely,  information 
on  the  boundary  layer  growth  constant,  C.  Additional  experi¬ 
mental  data  at  the  higher  Mach  numbers  was  presented  in  Refer¬ 
ences  (53)  and  (55),  but  at  a  time  so  that  they  could  not  be 
used  in  this  report.  These  new  experimental  data  should  be  in¬ 
corporated  into  Fig.  26  as  a  check  against  the  detailed  boundary 
layer  calculations  which  were  made  in  this  study  to  see  how  the 
new  data  fall  relative  to  the  low- stagnation  temperature  data 
taken  at  low  Mach  numbers.  The  above  studies  on  non-equilibrium 
flows,  isentropic  contours,  and  boundary  layer  correction  tech¬ 
niques  will  also  provide  information  on  the  constant,  C,  and  the 
entire  concept  of  the  approximate  boundary  layer  analysis. 
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FIG.  I  -  DRAG  COEFFICIENT  ON  A  CYLINDER  AT  M  =  5.92  UNDER 
RAREFIED  GAS  FLOW  CONDITIONS 
(FIG. 4  OF  REFERENCE  2) 
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FIG.2  -  SCHEMATIC  OF  LOW-DENSITY  WIND  TUNNEL 
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FIGURE  3  REYNOLDS  NUMBER  REQUIREMENTS  FOR  LOW 
DENSITY  WIND  TUNNELS 


STATIC  PRESSURE.p, MICRONS 


FIGURE  4-  VARIATION  OF  PRESSURE  WITH  MACH  NUMBER,  REYNOLDS 
NUMBER  AND  KNUDSEN  NUMBER  FOR  AIR,  To  =  7000R 
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REYNOLDS  NUMBER  PER  UNIT  LENGTH,  Re/L,  INCH"1 

FIGURE  5  -  VARIATION  OF  STATIC  PRESSURE  WITH  MACH  NUMBER,  REYNOLDS 
NUMBER  AND  KNUDSEN  NUMBER  FOR  AIR,  T0  -  I2000R 
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STAGNATION  PRESSURE,  pA  ,  ATM. 


ISENTRONIC 


(a)  ISENTROPIC  NOZZLE  CONTOUR 


(b)  FINAL  NOZZLE  CONTOUR  -  CORRECTED  FOR 
BOUNDARY  LAYER  GROWTH 

8*  -  BOUNDARY  LAYER  DISPLACEMENT  THICKNESS 
8  -  BOUNDARY  LAYER  THICKNESS  (99%  OF  FREE- 
STREAM  VELOCITY) 

Dse~  EXIT  DIAMETER  OF  ISENTROPIC  CONTOUR 
D^e  ACTUAL  EXIT  DIAMETER  OF  ISENTROPIC  CORE 
De  -  ACTUAL  EXIT  DIAMETE  OF  FINAL  CONTOUR 
D*  -  THROAT  DIAMETER 

FIG. 7“  NOMENCLATURE  FOR  APPROXIMATE 
BOUNDARY  LAYER  ANALYSIS 
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BOUNDARY  LAYER  THICKNESS  PARAMETER, 


FIGURE  8  BOUNDARY  LAYER  THICKNESS  PARAMETER 
FOR  A  SINE -ARCHVELOCITY  PROFILE 
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FIG.  9- COMPARISON  OF  LAMINAR  BOUNDARY  LAYER 

THICKNESS  FOR  SINE -ARCH  VELOCITY  PROFILE  AND 
FLAT  PLATE  SOLUTIONS 
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MACH  NUMBER, M 


FIG.  10- LAMINAR  BOUNDARY  LAYER  THICKNESS  CONSTANT 
FOR  COMPRESSIBLE  FLOW  OVER  A  FLAT  PLATE 
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MINIMUM  EXIT  DIAMETER  REYNOLDS  NUMBER 
FOR  ISENTROPIC  CONTOUR  ,  R«0## 


MACH  NUMBER, M 

FI6.II-  MINIMUM  EXIT  DIAMETER  REYNOLDS  NUMBER 
(ISENTROPIC  CONTOUR)  FOR  8/re  =  l 
AT  EXIT  OF  ACTUAL  NOZZLE 
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EXIT  DIAMETER  REYNOLDS  NUMBER 
(ISENTROPIC  CONTOUR)  Re0## 
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MACH  NUMBER, M 

FIG.  12-  DEPENDENCE  OF  EXIT  DIAMETER  REYNOLDS 

NUMBER  (ISENTROPIC  CONTOUR  )  ON  ISENTROPIC 
CORE  SIZE  AT  EXIT  OF  ACTUAL  NOZZLE 
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Re/L=  I03  INCH-1 
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MACH  NUMBER,  M 

FIGURE  13  ISENTROPIC  CONTOUR  EXIT  DIAMETER,  DISPLACEMENT,  THICKNESS  AND  ACTUAL  EXIT 
DIAMETER  FOR  A  LOW-DENSITY  WIND  TUNNEL  NOZZLE 
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ZERO  REACTION  RATE 
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AREA  RATIO,  A  /  A 

FIGURE  16  VARIATION  OF  COMPOSITION  FOR  THE  EXPANSION  OF  OXYGEN 
IN  A  15°  HYPERBOLIC  NOZZLE,  T0  =5900K,  p0  =  0.965  ATM. 
(FIG.  I  OF  REF.  35) 
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FIGURE  17  VARIATION  OF  COMPOSITION  FOR  THE  EXPANSION  OF  OXYGEN 
IN  A  15°  HYPERBOLIC  NOZZLE,  T0=5900  K,  p0*  9.4  ATM. 
(FIG.  I  OF  REF.  35) 
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FIG.  18-  VARIATION  OF  COMPOSITION  FOR  THE  EXPANSION  OF 
OXYGEN  IN  A  15°  HYPERBOLIC  NOZZLE, T0S 5900  K, 
D«=  82  ATM.  (FIG.  I  OF  REF.  35) 
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AREA  RATIO,  A/Aw 

FIG.  19- VARIATION  OF  STATIC  TEMPERATURE  FOR  THE 

EXPANSION  OF  OXYGEN  IN  A  15°  HYPERBOLIC  NOZZLE, 
T0  =  5900  K ,  p0  *  82  ATM.  (FIG.  2a  OF  REF.  35) 


APPENDIX  A 


STATIC  PRESSURE  REQUIREMENTS  FOR  LOW  DENSITY 
WIND  TUNNELS 


The  Knudsen  number  Is  defined  as: 


Kn  = 


(A-l) 


Multiplication  of  Eq.  (A-l)  by  several  ratios  yields: 


Kn  - 


V 


a 

o  V  L 
u 


(A- 2) 


Substituting  the  definitions  for  the  Mach  and  Reynolds  numbers 
into  Eq.  (A- 2)  yields: 


Kn  = 


M 

Re 


(A-3) 


For  the  pressures  and  temperatures  of  interest,  the  viscosity 
can  be  calculated  from  the  following  equation  (5) : 


where 


u  -33  o  v  X 


(A-4) 


(A- 5) 
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or 


v - 

TT  V 


(A-6) 

Substituting  Eqs.  (A-4)  and  (A-6)  into  Eq.  (A-3)  yields: 


Kn 


8  /  2  W 

5r  ii 


M_ 

Re 


A  speed  of  sound  parameter,  t  ,  is  defined  as: 


(A-7) 


2 

a _ p 


(A-8) 


The  compressibility  factor,  Z,  is  defined  as: 


W 

7  -  -E-  . 

0  T  R 


and 


zw  -  w 


Combining  Eqs.  (A-7)  through  (A-10)  yields: 


*-! •  fc 


For  flows  in  which  the  gas  composition  is  frozen: 


(A-9) 


(A-10) 


(A-ll) 


T  -  Y 


(A-12) 
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and  Eq.  (A-ll)  becomes: 


Kn 


1  ITT 

5  jJn 


M_ 

Re 


(A-13) 


The  Reynolds  number  per  unit  length  is  defined  as: 


Re  P  V  /a  i/\ 

- —  -  (A-14) 

L  u  v  ' 

Using  the  definition  of  the  Mach  number,  Eq.  (A-14)  becomes: 


Re 

L 


P  M  a 


(A-15) 


Substituting  Eos.  (A-8)  and  (A -9)  into  Ec.  (A-15)  yields: 


I  y 

Re  M  I  _o_  t_ 
L  * 


The  temperature  T  is  calculated  from: 


(A-16) 


=  1  +  m2 

T  2 


and  the  viscosity  can  be  calculated  from: 


(A-17) 


u  =  0.073225  x  10"5  T  ... 

201,6 

1  T 

Specific  units  are  involved  in  Ec.  (A-18),  namely, 

u  -  lbm/ft-sec 
T  -  °R 

1 16 


APPENDIX  B 


approximate  boundary  layer  analysis 

From  the  geometry  of  the  isentropic  contour  shown  in  Fig. 
7a,  the  following  expression  for  the  distance  between  the  throat 
and  the  exit  section  of  nozzle  can  be  derived: 


r* 

The  ratio  r  / r*  can  be  calculated  from  the  area  ratio  for  the 
se 

isentropic  contour,  or: 


The  area  ratio  is  known  from  the  specified  Mach  number  and  value 
of  v  for  a  particular  case.  Equation  (B-l)  can  be  rewritten  as: 


L 


t 


r 


se 


(B-3) 


It  is  assumed  that  the  boundary  layer  thickness  at  the  exit 
of  the  nozzle  can  be  calculated  from  the  following  equation: 


(B-4) 
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where  C  is  obtained  from  Fig.  26.  Equation  (B-4)  can  be  re¬ 
written  as: 


c  AT  V°Te 


(B-5) 


where  the  exit  diameter  Reynolds  number  for  the  isen tropic 


contour.  Re 


Dse 


is  defined  as: 


Re 


p  V  D 
se  se  se 


Dse 


u 


se 


Combining  Eqs.  (B-3)  and  (B-5)  yields 

6  F  C 


J  COC  ~  tjl 


se 


4 


Re 


Dse 


(B-6) 


(B-7) 


The  actual  exit  radius  of  the  nozzle,  rg,  is  determined  by  the 
equation: 


r  -  r  +6*  (B-8) 

e  se  '  7 

Since  the  boundary  layer  thickness  and  the  displacement  thick¬ 
ness  are  related  by  the  following  equation: 


6*  -  6(1  "  B)  (B-9) 

Equation  (B-8)  can  be  written: 

r  -  r  +  (1  -  B)  6  (B-10) 

e  se 
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Thus  an  expression  for  6/r  is: 

e 


r  r  (1  -  B)6 
e  se  v  J 


(rge/6)  X  (1  -  B) 


or  substituting  Eq.  (B-7)  into  Eq.  (B-ll)  yields: 


(B-ll) 


r 


(B-12) 


+  (1  -  B) 


The  following  expression  for  Re^se  can  be  obtained  from 


Eq.  (B-12): 


_  2  2  e 

Re_  =  2  B  C  cot  -r 
Dse  2 


[x  •  fifcnTT  '  x]  +  1 J 


(B-13) 

Equation  (B-5)  can  be  combined  with  Eq.  (B-3)  to  give: 


/— 1 

r-rr 

ft  -  c 

Jr 

J  cot  -r  /  1 

1  se 

V  2  ^ 

1 


Further  rearrangement  of  Eq.  (B-14)  yields: 


(B-14) 
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or 


/ ?c  r~  e 1  r  i  * 
—  lcot  2  J1  -  firm 


ReDse 
<ReDse/Dse> 


(B-16) 


However,  Eq.  (B-13)  is  an  expression  for  ReDge  and: 


t221  -  r  (B-17) 

se 

Using  this  information  in  Eq.  (B-16)  yields  the  final  expression 
for  6 : 

“2  ■  cot  f  •  [  1  [(77T)  -  ]  *  '} 

5  ■  - ^7l - 

(B-18) 
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FIGURE  20  VARIATION  OF  STATIC  TEMPERATURE  FOR  THE  EXPANSION  OF  AIR 
IN  A  HYPERSONIC  NOZZLE,  T.-6000K,  p.-IOOATM. 

(FIGURE  29e  OF  REF.  20) 
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AREA  RATIO,  A/A# 

FIGURE  22  VARIATION  OF  VELOCITY  FOR  THE  EXPANSION  OF  AIR  IN  A  HYPER¬ 
SONIC  NOZZLE  To  =  6000K,  p0=IOOATM. 

(FIGURE  29  d  OF  REF.  20) 


AREA  RATIO,  A/Aw 

FIGURE  23  VARIATION  OF  MACH  NUMBER  FOR  THE  EXPANSION  OF  AIR  IN  A 
HYPERSONIC  NOZZLE,  To*6000K,  p0=IOOATM. 

(FIGURE  29c  OF  REF.  20) 
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AREA  RATIO,  a/a* 

FIGURE  24  VARIATION  OF  STATIC  PRESSURE  FOR  THE  EXPANSION  OF  AIR  IN  A 
HYPERSONIC  NOZZLE,  T0»6000K,  p0»  IOO  ATM. 

(FIGURE  29b  OF  REF.  20) 
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FIGURE  25  VARIATION  OF  DENSITY  FOR  THE  EXPANSION  OF  AIR  IN  A  HYPERSONIC 
NOZZLE ,  T0  =6000K  ,  p0=  IOO  ATM. 

(FIGURE  29a  OF  REF.  20) 


r* 

1.4 

r »  1.2 

EXPERIMENTAL 

CODE 

REF 

A 

SECOND  NOZZLE 

O 

(14) 

DESIGN  THIS  REPORT 

ft 

(15) 

© 

(16) 

3 

(17) 

A 

THIRD  NOZZLE 

c 

(16) 

DESIGN  THIS  REPORT 

O 

(19) 

THEORETICAL 
O  FIRST  NOZZLE 


DESIGN  THIS  REPORT 


MACH  NUMBER, M 


FIGURE  26  RECOMMENDED  VALUES  FOR  THE  BOUNDARY 
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MASS  FLOW  RATE,  w/A  ,  LBM/SEC-IN 


|  10  I02  I03 

STAGNATION  PRESSURE,  p0  ,  ATMOSPHERES 

FIG.  27-  MASS  FLOW  RATE  REQUIREMENTS  FOR  LOW- DENSITY 
WIND  TUNNELS  USING  AIR 
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FIG. 28  HEATER  REQUIREMENTS  FOR  LOW- DENSITY  WIND 
TUNNELS  USING  AIR 


FIG. 29.  VOLUME  FLOW  RATE  REQUIREMENTS 
FOR  LOW- DENSITY  WIND 
TUNNELS  USING  AIR 
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